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This  AFOSR  YIP  Program  sought  to  investigate  the 
fabrication  and  properties  of  3D  “morphing”  carbon  nanotube 
(CNT)  microstructures.  This  concept  was  inspired  by  the  design 
and  responsive  behavior  of  plant  tissues  that  swell  or  contract  in 
response  to  a  stimulus,  and  it  was  proposed  that  morphing  CNT 
microstructures  would  operate  by  swelling  of  an  active  material  in 
a  direction  locally  perpendicular  to  the  CNTs  (Fig.  1).  This  shape 
change,  in  combination  with  the  mechanical  and  electrical 
characteristics  of  the  CNTs,  could  enable  incorporation  of  intrinsic 
actuation  and  sensing  behaviors  in  design  of  new  structural 
materials  and  active  surfaces. 

This  final  performance  report,  presented  in  four  self-contained 
chapters,  describes  the  major  accomplishments  of  the  AFOSR 
program.  These  include  advances  in  the  fabrication  and 
mechanics  of  CNT  microstructures,  and  demonstrations  of  novel 
composite  materials  and  active  materials  based  on  composites  of 
organized  CNTs,  polymers,  and  ceramics.  The  chapters  are: 

1:  Corrugated  paraffin  nanocomposite  films  as  large 

stroke  thermal  actuators  and  self-activating  thermal 

interfaces.  We  demonstrate  that  compact,  high- 
performance  paraffin  actuators  can  be  made  by  confining 
paraffin  within  vertically  aligned  carbon  nanotube  films. 

This  large-stroke  vertical  actuation  is  enabled  by  strong 
capillary  interaction  between  paraffin  and  CNTs  and  by 
engineering  the  CNT  morphology  by  pre -compression  before  capillary-driven  infiltration  of  the 
molten  paraffin.  The  maximum  actuation  strain  (-0.02-0.2)  approaches  that  of  natural  muscle  and 
exceeds  that  of  shape  memory  alloys.  Further,  we  show  how  a  CNT -paraffin  film  can  serve  as  a 
self-activating  thermal  interface  that  closes  a  gap  when  it  is  heated.  These  new  CNT -paraffin 
film  actuators  could  be  scaled  to  large-area  production  and  lamination,  and  are  attractive  for 
practical  use  in  miniature  systems  due  to  their  self-contained  design. 

2:  Replica  molding  of  liquid  crystal  network  (LCN)  polymer  microstructures  for  active 

surfaces.  We  present  a  replica  molding  (RM)  fabrication  technique  to  obtain  high  fidelity 
surface  bound  glassy  LCN  microstructure  actuators  over  large  areas.  The  removal  of  oxygen 
during  curing  was  critical  to  achieve  high-fidelity  replicas  with  smooth  surface  texture,  while  the 
necessary  LCN  alignment  for  actuation  was  obtained  via  curing  in  a  magnetic  field.  Cast 
microstructures  had  lateral  dimensions  of  a  few  microns  and  aspect  ratios  (AR)  up  to  15: 1 .  A 
gradient  based  edge  detection  algorithm  was  used  to  quantify  the  LCN  microstructure  actuation, 
verifying  anisotropy  of  the  network. 

3:  Strain-engineered  manufacturing  of  freeform  carbon  nanotube  microstruetures.  Here, 

we  demonstrate  a  new  technique  for  scalable  manufacturing  of  freeform  microstructures,  via 
strain-engineered  growth  of  aligned  carbon  nanotubes  (CNTs).  Offset  patterning  of  the  CNT 
growth  catalyst  is  used  to  locally  modulate  the  CNT  growth  rate.  This  causes  the  CNTs  to 
collectively  bend  during  growth,  with  exceptional  uniformity  over  large  areas.  The  final  shape  of 
the  curved  CNT  microstructures  can  be  designed  via  finite  element  modeling,  and  compound 
catalyst  shapes  produce  microstructures  with  multi -directional  curvature  and  unusual  self- 
organized  patterns.  Conformal  coating  of  the  CNTs  enables  tuning  of  the  mechanical  properties 
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Fig.  1.  Morphing  CNT  microstructure 
concept:  (a)  lateral  swelling  of  aligned 
CNT  network;  (b)  radial  expanding 
actuator;  (c)  lateral  bending  actuator; 

(d)  vertical  twisting  actuator. 
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independently  from  the  microstructure  geometry,  representing  a  versatile  principle  for  design  and 
manufacturing  of  complex  microstructured  surfaces. 

4:  Scaling  the  stiffness,  strength,  and  toughness  of  ceramic-coated  nanotube  foams  into  the 
structural  regime.  We  demonstrate  how  to  use  aligned  carbon  nanotube  (CNT)  growth  followed 
by  atomic  layer  deposition  (ALD)  to  manufacture  nanocomposite  foams  with  precisely 
controllable  mechanical  properties.  By  starting  with  a  low  density  CNT  forest  and  varying  the 
ALD  coating  thickness,  we  realize  predictable  ~1000-fold  control  of  Young’s  modulus  (14  MPa 
to  20  GPa,  where  E~p2.8),  ultimate  compressive  strength  (0.8  MPa  to  0.16  GPa),  and  energy 
absorption  (0.4  to  400  J/cm3).  Owing  to  the  continuous,  long  CNTs  within  the  ceramic 
nanocomposite,  the  compressive  strength  and  toughness  of  the  foam  are  10-fold  greater  than 
commercially  available  aluminum  foam  over  the  same  density  range.  Moreover,  the  compressive 
stiffness  and  strength  equal  that  of  compact  bone  at  10%  lower  density.  Along  with  emerging 
technologies  for  scalable  patterning  and  roll-to-roll  manufacturing  and  lamination  of  CNT  films, 
coated  CNT  foams  may  enable  new  materials  for  mechanical  and  multifunctional  applications 
such  as  catalysis,  filtration,  and  thermal  protection. 

Future  work  may  focus  on  achieving  extreme  stimuli -responsive  deformations  using  CNT-polymer 
and  CNT-LCN  systems,  and  enabling  multifunctional  surfaces  combining  the  strain-dependent  electrical 
properties  of  CNTs  with  their  active  response. 
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Chapter  1:  Corrugated  paraffin  nanocomposite  films  as  large  stroke  thermal  actuators  and 
self-activating  thermal  interfaces 

Contents  of  this  chapter  have  been  submitted  for  publication  as:  D.  Copic  and  A.J.  Hart, 
“Corrugated  paraffin  nanocomposite  films  as  large  stroke  thermal  actuators  and  self-activating 
thermal  interfaces.  ” 

Introduction 

The  development  of  active  materials  whose  dimensions  can  be  changed  using  externally 
applied  electrical,  optical,  or  thermal  stimuli,  is  essential  to  advances  in  actuators  for  use  in 
miniaturized  grippers  and  pumps,  probes  for  scanning  probe  microscopy,  and  ultrasonic 
transducers.  Specifically,  paraffin  wax  has  found  longstanding  use  in  linear  actuators  with 
centimeter-scale  and  larger  dimensions,  which  are  widely  used  in  automotive  and  aerospace 
systems.  This  is  because  paraffin  wax  exhibits  a  large  and  repeatable  volume  change  (10-40%) 
and  stress  output  (10  -  100  MPa)1-2  upon  melting,  resulting  in  a  high  mechanical  energy  density 
(107  J  m"3)2  (Figure  SI).  In  comparison  to  paraffin,  thin  film  shape  memory  alloy  (SMA) 
actuators  exhibit  actuations  strains  of  7x10  3  -  7x10  2  and  stresses  of  100  -  700  MPa,  while  thin 
film  piezoelectric  (PZT)  actuators  exhibit  actuation  strains  of  5xl0"6  -  2x1 0“4  and  stresses  of  1  - 
9  MPa3. 

Paraffin  has  also  been  incorporated  in  micro-scale  systems  including  on-chip  pumps  4,  valves 
5_8,  thin-film  cantilevers  5-9,  and  linear  actuators10.  These  designs  and  the  related  fabrication 
methods  are  complicated  by  the  need  to  confine  the  paraffin  during  melting,  which  is  needed  to 
convert  the  expansion  of  the  paraffin  wax  into  a  directed  motion.  As  a  result,  the  overall  system 
performance  of  miniature  paraffin  actuators  is  often  compromised  relative  to  the  intrinsic 
actuation  capability  of  paraffin,  and  these  actuators  require  a  large  footprint  compared  to  the  size 
of  active  material  used. 

Therefore,  owing  to  its  intrinsic  performance  and  simplicity  of  operation,  it  would  be  useful 
to  develop  more  scalable  means  to  harness  the  actuation  capability  of  paraffin  wax  in  composite 
materials.  Recently,  Lima  et  al.  demonstrated  micro-fiber  actuators  comprised  of  carbon 
nanotube  (CNT)  yarns  infiltrated  with  paraffin  wax  11,12.  In  these  actuators,  the  wax  melts  inside 
the  yarn  upon  heating,  and  remains  confined  within  the  yarn  due  to  its  strong  wetting  of  the 
CNTs  and  high  surface  tension.  As  the  wax  expands,  it  swells  the  CNT  yarn,  causing  the 
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composite  to  contract  its  length  due  to  the  mechanical  anisotropy  introduced  by  the  twist  of  the 
CNTs  about  the  yarn  axis.  These  yarns  exhibited  up  to  8%  contraction  upon  heating  11  and  were 
able  to  lift  17,700  times  their  own  weight  with  a  3%  contraction. 

Here,  we  show  how  large  stroke  thin  film  nanocomposite  actuators  can  be  made  by  confining 
paraffin  within  vertically  aligned  CNT  films  (CNT  “forests”)  and  how  the  actuation  stroke  of  the 
composite  can  be  engineered  by  corrugating  the  forests.  Via  this  strategy,  we  demonstrate  the 
fabrication  of  large-stroke  CNT-paraffin  composite  film  actuators,  and  show  that  the  actuation 
performance  is  governed  by  two  distinct  mechanistic  regimes  of  operation:  (1)  low  strain 
expansion  below  the  wax  melting  point,  wherein  the  composite  actuator  has  high  stiffness  and 
(2)  high  strain  expansion  during  and  beyond  the  melting  transition,  wherein  the  actuator  has 
much  lower  stiffness.  Last,  we  show  how  the  high  strain  regime  can  be  exploited  to  build  a  self¬ 
activating  thermal  interface,  where  the  vertical  extension  of  the  nanocomposite  film  causes  it  to 
contact  a  second  substrate,  resulting  in  enhanced  heat  dissipation  at  elevated  temperature. 

Methods 

CNT  growth:  CNT  forests  were  grown  using  chemical  vapor  deposition  (CVD)  from  a  thin  film 
catalyst  (1  nm  Fe  supported  on  10  nm  AI2O3),  which  was  sputtered  on  a  (100)  Si  wafer  coated 
with  300  nm  thermally  grown  SiCh.  The  wafer  was  cut  into  4  mm  by  8  mm  chips  using  a  dicing 
saw.  Two  chips  were  placed  side-by-side  inside  a  horizontal  tube  furnace  (Scientific  Mini-Mite) 
at  atmospheric  pressure.,  which  was  then  flushed  with  helium,  heated  in  a  400:100  seem  H2:He 
gas  mixture  for  10  minutes  to  775°C,  annealed,  held  under  these  conditions  for  a  further  10 
minutes,  and  exposed  to  a  100:400:100  seem  H2:He:C2H4  gas  mixture  for  20-22  minutes  for 
CNT  growth.  Rapid  cooling  of  the  CNT  forests  was  achieved  by  shifting  the  quartz  tube 
downstream  until  the  CNT  forests  were  outside  the  furnace. 

Paraffin  wax  infiltration:  CNT  forests  were  inverted  such  that  their  top  made  contact  with  a  2cm 
by  2cm  silicon  chip,  which  was  heated  to  100°C  by  a  hotplate.  The  forests  were  then  compressed 
by  10-50%  of  its  initial  height  using  a  vertical  precision  actuator.  While  maintaining  the  position 
and  temperature,  molten  paraffin  wax  was  added  to  the  silicon  chip  and  promptly  wicked  into  the 
CNT  forest.  Constant  load  was  applied  during  sample  cooling  and  the  sample  was  manually 
separated  from  the  silicon  chip.  The  sample  was  then  placed  on  a  spin  coater  and  exposed  it  to 
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75mW/cm2  of  ultra  violet  (UV)  light  (Dymax  Model  2000)  for  35s  prior  to  spinning.  Excess 
paraffin  was  then  spun  away  at  3000  rpm  for  30  seconds  while  the  UV  light  was  maintained. 
Displacement  and  load  measurements:  All  forest  indentation,  displacement,  and  force 
measurements  under  load  were  collected  using  a  custom-built  single  axis  precision  actuator  with 
both  PID  position  and  load  control.  The  custom  precision  actuator  was  driven  by  a  voice  coil, 
while  the  actuator  motion  was  measured  using  an  optical  encoder  and  the  load  was  measured 
using  a  S-beam  load  cell  (Futek).  The  actuator  displacement  under  no  load  was  obtained  by  an 
edge-tracking  algorithm,  which  was  used  on  a  series  of  high-resolution  images  taken  by  a 
telecentric  lens  of  the  actuator  side  view  backlit  by  a  collimated  white  light  source. 

Thermal  conductivity  measurement:  CNT-paraffin  actuators  were  placed  on  a  PI  controlled 
hotplate  (Torrey  Pines  Scientific  HS40)  instrumented  with  a  heat  flux  sensor  and  thermocouple 
(Omega  HFS-4),  while  a  similarly  instrumented  heat  sink  was  separated  from  the  actuator  by 
100pm  using  a  micrometer  stage.  The  temperature  of  the  hotplate  was  slowly  ramped  from  room 
temperature  to  80°C  while  recording  the  heat  fluxes  and  temperatures.  Displacement 
measurements  were  obtained  using  the  edge  tracking  algorithm  and  setup  described  above.  The 
thermal  resistivity  was  calculated  using  the  measured  data  and  a  simple  thermal  network 
resistance  model,  which  accounted  for  the  convective  heat  transfer  between  the  heat  source  and 
sink,  the  thermal  contact  resistance,  and  the  resistance  of  the  silicon  substrate  and  air  gap. 

Results  and  Discussion 

CNT-paraffin  nanocomposite  films  were  fabricated  by  capillary  driven  infiltration  of  molten 
paraffin  wax  into  vertically  aligned  CNT  forests,  as  shown  in  Figure  1.  First,  CNT  forests  were 
grown  on  silicon  wafers  by  atmospheric  pressure  thermal  chemical  vapor  deposition  (CVD) 
using  a  supported  catalyst  (1/10  nm  Fe/AFCF)  deposited  by  sputtering  (see  Methods)13,14.  After 
the  allotted  growth  time,  which  determines  the  CNT  forest  height,  the  CVD  system  was  cooled 
rapidly  while  maintaining  the  hydrocarbon  flow;  this  resulted  in  strong  CNT-substrate  adhesion, 
which  prevented  delamination  of  the  forest  during  paraffin  expansion. 

Next,  the  CNT  forest  was  compressed  vertically  using  a  precision  actuator,  while  resting  on  a 
hot  plate.  Upon  compression  the  CNTs  collectively  folded,  starting  at  the  base  of  the  forest, 
resulting  in  the  corrugated  morphology  shown  in  Figure  ID15  18.  The  corrugated  morphology  is 
visually  reminiscent  of  a  macroscopic  accordion  bellows.  Unless  otherwise  noted  below,  the 
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CNT  forest  was  compressed  to  50%  of  its  as-grown  height.  Next,  the  corrugated  CNT  forest  was 
heated  to  100°C  and  molten  paraffin  wax  was  poured  onto  to  the  sample  from  the  side,  while 
maintaining  the  compressive  load.  The  wax  rapidly  wetted  the  CNTs  and  wicked  into  the  forest. 
The  hot  plate  was  then  cooled  (~4°C/min)  to  35 °C  while  keeping  the  sample  under  constant  load. 
Last,  the  wax-infiltrated  CNT  forest  was  removed  from  the  hot  plate,  and  placed  on  a  spin- 
coater.  The  sample  was  heated  using  a  UV  lamp,  which  caused  the  wax  to  melt,  and  spinning 
was  used  to  remove  excess  wax  from  the  substrate  while  retaining  wax  within  the  CNT  forest 
due  to  capillary  action. 

The  corrugated  forest  morphology  is  essential  to  achieving  large  actuation  stroke  upon 
paraffin  phase  change.  Maintaining  compressive  load  during  paraffin  infiltration  and  cooling 
allows  the  corrugated  morphology  to  be  maintained  in  the  nanocomposite  film.  If  the  CNT 
forest  is  unloaded  before  paraffin  infiltration,  the  wave-like  corrugation  is  only  partially 
recovered  because  the  elastic  energy  stored  in  the  deformed  CNTs  is  counteracted  by  the  van  der 
Walls  adhesion  between  CNT  contacts  created  during  forest  compression.  If  the  forest  is 
infiltrated  with  paraffin  without  maintaining  the  compressive  load,  the  forest  expands  to 
approximately  80%  of  its  initial  height  during  cooling  and  develops  cracks  in  its  top  surface 
during  cooling  (Figure  S2).  These  cracks,  which  extend  through  to  the  substrate,  are  caused  by 
local  capillary  aggregation  of  the  CNTs,  which  is  driven  by  the  shrinkage  of  the  wax  upon 
cooling.  The  pressure  from  capillary  aggregation  exceeds  the  lateral  strength  of  the  forest  and 
causes  these  cracks.  Upon  thermal  cycling,  such  cracked  composite  films  close  and  open  their 
cracks,  instead  of  expanding  and  contracting  vertically,  and  therefore  do  not  operate  as  intended. 

Upon  heating,  the  crack-free  nanocomposite  film  expands  vertically  and  the  wavelength  of 
the  corrugations  increases,  visually  resembling  inflation  of  an  accordion  bellows.  A  side  view 
optical  image  at  25 °C  and  a  mirrored  image  of  the  same  film  heated  to  150°C  are  shown  in 
Figure  2A.  Measurements  of  the  thermally  induced  strain  with  no  load  applied  were  made  via 
edge  tracking  of  optical  images  taken  using  a  telecentric  lens  backlit  by  a  collimated  white  light 
source.  This  actuator  extended  vertically  by  -17%  (perpendicular  to  the  substrate)  at  150°C,  as 
shown  in  Video  S3.  Strain  versus  temperature  data  from  three  consecutive  heating  and  cooling 
cycles  (Figure  2B)  shows  that  the  vertical  actuation  of  the  nanocomposite  is  repeatable. 

With  increasing  temperature  below  the  melting  point,  the  CNT-paraffin  film  expands  slowly, 
then  it  expands  rapidly  through  the  melting  transition  (53-57°C),  and  after  that  it  expands  in  an 
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approximately  linear  relationship  with  temperature.  Thus,  the  strain-temperature  curve  of  a 
characteristic  CNT-paraffin  actuator  has  three  distinct  regions:  (1)  low  thermal  expansion  below 
the  melting  point,  (2)  rapid  expansion  during  melting,  and  (3)  further  expansion  driven  by  the 
increased  pressure  of  the  confined  molten  wax.  In  the  second  and  third  regions,  the  shape  of  the 
curve  for  the  CNT-paraffin  nanocomposite  is  similar  to  that  for  unconstrained  paraffin  wax  of 
similar  molecular  weight  19 ,  but  exhibits  a  lower  volume  change.  Through  the  complete  thermal 
cycle,  the  CNT-paraffin  films  undergo  a  volume  change  of  approximately  20%,  compared  to 
approximately  40%  expected  for  unconstrained  neat  paraffin  with  similar  molecular  weight.  The 
observed  volume  change  is  reduced  because  the  paraffin  must  work  against  the  stiffness  of  the 
CNT  forest  to  expand  the  composite.  Comparable  mechanisms  limit  the  strain  of  macroscale 
paraffin  actuators,  yet  to  a  smaller  degree,  wherein  the  paraffin  is  held  within  a  container.  For 
further  discussion,  we  denote  the  low  thermal  expansion  region  below  the  melting  point  as 
“regime  1”  and  the  subsequent  behavior  of  rapid  expansion  upon  melting  and  linear  expansion 
beyond  melting  as  “regime  2”. 

To  explore  the  relationship  between  applied  force  and  actuation  strain,  we  applied  constant 
loads  during  heating  while  monitoring  the  vertical  displacement  of  the  film,  using  a  feedback- 
controlled  loading  system.  Consecutive  thermal  strain  curves  under  constant  stress  levels  varying 
by  factors  of  ten,  from  0.01  kPa  to  10  kPa,  are  shown  in  Figure  2C.  In  regime  1  (below  the  wax 
melting  point),  all  four  strain  curves  overlap  closely  and  reach  a  maximum  thermal  strain  of  3% 
at  50°C.  For  all  but  the  lOkPa  curve,  rapid  expansion  upon  melting  takes  place.  As  expected,  an 
increase  in  applied  load  leads  to  smaller  actuation  strains.  At  lOkPa  load,  the  CNT-paraffin 
nanocomposite  collapses  after  melting,  suggesting  that  the  blocking  stress  of  the  CNT-paraffin 
actuator  in  regime  2  (above  the  wax  melting  point)  is  between  1- lOkPa.  In  the  data  shown,  the 
strains  reached  under  load  surpass  those  reached  in  the  no  load  condition,  which  is  likely  due  to 
run-to-run  variation  in  CNT  forest  density20. 

Once  the  wax  is  melted  inside  the  CNT  forest  (>60°C),  the  vertical  stiffness  of  the  film 
decreases  significantly,  as  shown  in  Figure  3.  This  signifies  the  onset  of  regime  2.  At  room 
temperature  the  modulus  of  the  actuator  is  18.6  MPa,  which  is  approximately  15-fold  greater 
than  a  corrugated  (crushed)  CNT  forest  that  has  not  been  infiltrated  with  paraffin.  The  actuator 
stiffness  unexpectedly  increases  to  46.7  MPa  at  40°C.  The  actuator  stiffness  then  drops 


5 


significantly,  to  4.5MPa,  as  the  wax  melts  at  60°C,  while  at  100°C  the  stiffness  of  the  actuator 
drops  to  1.2  MPa,  which  is  less  than  the  stiffness  of  the  corrugated  CNTs  without  paraffin. 

In  regime  1,  actuation  of  the  CNT-paraffin  film  is  powered  by  the  release  of  the  deformation 
energy  (stored  in  the  CNTs  during  wax  freezing),  which  does  work  against  the  viscous  forces  in 
the  wax  and  lifts  the  applied  load  upwards  by  up  to  3%  strain  at  lOkPa.  Because  the  compressive 
load  used  to  crush  the  forest  is  maintained  during  cooling  of  the  freshly  infiltrated  was,  and 
because  wax  undergoes  shrinkage  during  cooling,  additional  elastic  energy  is  stored  in  the  CNT 
network  during  freezing  of  the  wax.  As  the  wax  softens  before  melting,  it  allows  the  corrugated 
microstructure  to  relax,  releasing  the  stored  deformation  energy.  When  the  wax  melts,  the 
wavelength  of  the  corrugated  CNT  morphology  increases,  signaling  the  onset  of  regime  2.  The 
behavior  of  the  composite  in  regime  1,  below  the  wax  melting  point,  is  analogous  to  shape 
memory  materials,  where  potential  energy  is  first  stored  during  cooling  and  is  later  released 
during  heating,  resulting  in  a  shape  change.  However,  unlike  most  shape  memory  materials, 
upon  cooling  the  wax  pulls  the  CNTs  back  to  their  original  shape,  without  requiring  an  external 
restoring  force.  Therefore,  the  CNT-paraffin  system  is  a  fully  reversible  shape  memory 
composite,  and  to  our  knowledge  this  is  the  first  time  that  such  has  been  achieved  using  paraffin 
wax  as  the  active  material. 

In  regime  2,  the  wax  has  undergone  melting  and  expands  up  to  20%  more  in  total,  doing 
work  to  extend  the  corrugated  CNT  forest.  The  equilibrium  position  of  the  actuator,  at  elevated 
temperature,  is  reached  when  pressure  of  the  molten  wax  inside  the  forest  equals  the  pressure 
required  to  incrementally  extend  the  forest  in  the  vertical  direction  in  addition  to  any  externally 
applied  load.  If  the  pressure  exerted  by  the  paraffin  exceeds  the  capillary  pressure,  the  molten 
wax  is  expected  to  flow  out  of  the  forest.  For  example,  loading  an  actuator  to  lOOkPa  at  150°C 
results  in  the  escape  of  wax  and  the  film  becomes  compressed  to  less  than  10%  of  its  initial 
height. 

The  actuation  properties  of  CNT-paraffin  films  are  strongly  dependent  on  the  compression 
strain  applied  to  the  CNT  forest  before  and  during  paraffin  infiltration;  we  found  that  samples 
compressed  by  approximately  50%  of  their  initial  height  resulted  in  the  actuators  with  the  largest 
vertical  stroke  (Figure  S4).  When  compressed  by  50%  of  its  as-grown  thickness,  a  CNT  forest 
was  observed  to  have  a  corrugated  conformation  for  the  lower  45-50%  of  its  compressed 
thickness  (as  viewed  from  the  side).  Further,  the  density  and  alignment  of  the  CNTs  within  the 
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forest  are  expected  to  influence  the  capillary  pressure  confining  the  wax  within  the  composite, 
therefore  changing  the  actuator  stiffness  and  strain  output.  A  greater  CNT  density  would  result 
in  a  lower  average  spacing  between  individual  CNTs,  resulting  in  a  greater  capillary  pressure  of 
the  paraffin  and  therefore  a  greater  blocking  stress  of  the  actuator.  Additionally,  we  expect  a 
greater  CNT  forest  density  will  correspond  to  greater  actuator  stiffness  at  the  expense  of  less 
maximum  strain.  For  regime  2,  this  tradeoff  is  due  to  the  energy  density  of  the  process  powering 
the  shape  change,  in  this  case  the  phase  transition  and  thermal  expansion  of  paraffin,  remaining 
largely  unchanged  by  the  density,  morphology,  or  increased  stiffness  of  the  CNT  forests. 

While  the  CNT-paraffin  actuators  exert  a  relatively  low  stress  after  melting  (clOkPa),  the 
maximum  actuation  strain  (-0.02-0.2)  approaches  that  of  natural  muscle  (-0.3-0. 7)  and  exceeds 
that  of  shape  memory  alloys  (-0.07).  Further  study  is  needed  to  explore  the  mechanism  and 
limits  of  actuation,  yet  it  is  clear  that  the  strain  and  stress  will  be  related  to  the  density  of  the 
CNTs,  their  initial  deformation,  and  the  intrinsic  active  behavior  of  the  paraffin  wax.  Moreover, 
although  the  blocking  stress  for  regime  1  is  not  known,  the  present  CNT-paraffin  film  actuator  is 
able  to  lift  at  least  2,000  times  its  weight  at  100°C.  A  comparison  of  the  CNT-paraffin  film 
performance  to  traditional  actuators  and  active  materials  can  be  found  in  Figure  S6. 

Last,  we  explore  one  potential  application  of  the  CNT-paraffin  film  actuator,  a  passive 
thermal  switch21,  which  enables  engineered  heat  dissipation  via  self-actuated  contact  with  a 
second  surface.  We  present  a  preliminary  demonstration  of  this  by  measuring  the  thermal 
conductivity  of  the  CNT-paraffin  actuator  as  it  expands  upward  to  contact  a  second  surface.  We 
constructed  a  setup  with  thin  film  heat  flux  sensors  on  parallel  surfaces  (Figure  4A),  and  the 
measured  heat  flux  through  the  composite  versus  temperature  as  it  closes  the  gap  between  the 
surfaces  (Figure  4B).  The  thermal  conductivity  of  the  CNT-paraffin  composite  was  calculated 
from  a  combined  thermal  resistance  model  using  the  measured  heat  flux  into  the  heat  sink  and 
the  temperature  difference  between  the  hotplate  and  the  heat  sink  (see  SI  for  model  description). 

When  the  system  is  heated  from  below  the  CNT-paraffin  film,  the  thermal  resistance  is 
initially  large  while  most  of  the  heat  is  absorbed  by  the  paraffin  wax.  As  the  film  expands,  the 
thermal  resistance  of  the  actuator  and  the  gap  decreases  rapidly,  the  air  gap  is  reduced,  the 
conductivity  of  the  paraffin  increases  22,  and  ultimately  the  actuator  contacts  the  heat  sink.  As 
the  gap  closes  and  the  composite  contacts  the  upper  heat  sink,  the  thermal  conductivity  increases 
to  approximately  0.4  W/m-K  at  70°C  (Figure  S5).  This  is  approximately  twice  that  of  neat 
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paraffin  (0.21-0.24  W/m-K).  The  greater  thermal  conductivity  than  neat  paraffin  is  due  to  the 
presence  of  aligned  CNTs  within  the  composite  (~1  wt%).  The  final  value  is  comparable  to 
commercially  available  polymer-ceramic  composite  thermal  transfer  tapes23,  and  is  close  to  that 
of  passive  composites  made  by  polymer  infiltration  to  CNT  forests24,25.  In  principle,  the 
temperature  at  which  the  actuator  engages  the  heat  sink  could  be  controlled  via  the  initial  gap  or 
the  physical  characteristics  of  the  CNT-paraffin  film  (pre-strain,  paraffin  composition,  CNT 
density,  etc.). 

Conclusion 

In  conclusion,  we  have  shown  that  large-stroke  thermal  actuators  can  be  made  by  infiltration 
of  corrugated  CNT  forests  with  paraffin  wax.  Two  mechanistic  regimes  of  the  CNT-paraffin 
actuator  were  identified:  (1)  a  high  stiffness,  low  strain  regime  below  the  melting  point,  which  is 
analogous  to  the  performance  of  shape  memory  materials  and  (2)  a  low  stiffness,  high  strain 
regime  above  the  melting  point.  In  regime  1  the  actuator  was  capable  of  lifting  2,000  times  is 
own  weight  with  an  actuation  strain  of  3%  and  a  stress  of  lOkPa  at  40°C,  while  in  regime  2  the 
actuator  was  shown  to  reach  an  actuation  strain  of  20%  and  a  stress  of  lkPa  at  175°C.  This 
large-stroke  vertical  actuation  is  enabled  by  strong  capillary  interaction  between  paraffin  and 
CNTs,  and  engineering  of  the  vertical  compliance  by  corrugation  of  the  CNT  forest.  CNT- 
paraffin  nanocomposites  have  been  shown  to  act  as  switchable  thermal  interfaces  with  a  thermal 
conductivity  up  to  0.4  W/m-K  at  70°C  with  a  CNT  content  of  1  wt%.  Because  the  CNT-paraffin 
film  actuators  do  not  require  external  confinement  to  restrain  the  wax,  the  materials  are  amenable 
to  integration  with  laminates  and  miniaturization  for  applications  including  thermal  switching 
and  flow  control. 
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Figure  1.  Fabrication  of  corrugated  CNT-paraffin  nanocomposite  films:  (a-c)  Key  process  steps; 
SEM  images  of  (d)  compressed  CNT  forest  with  corrugations  and  (e)  CNT-paraffin  composite 
film,  with  (f)  close-up  of  corrugated  microstructure. 
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Figure  2.  Actuation  of  CNT-paraffin  composite  films:  (a)  Image  of  the  actuator  side  view  at 
25  °C  and  a  mirrored  image  of  the  same  view  at  150°C  with  1mm  scale  bar.  (b)  Three 
consecutive  cycles  of  a  CNT-paraffin  actuator  under  no  load,  (c)  Displacement  and  strain  versus 
temperature  under  varying  values  of  constant  applied  stress  during  heating. 


13 


a 


60 


20  40  60  80  100 


Temperature  (”C) 


Figure  3.  (a)  Temperature  dependent  compressive  stiffness  of  CNT-paraffin  films  and  (b) 
exemplary  stress-strain  curves  used  to  determine  stiffness. 
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Figure  4.  Demonstration  of  CNT-paraffin  actuator  as  a  passive  thermal  switch:  (a)  Simplified 
setup  of  measurement  configuration;  (b)  measured  heat  flux  through  the  actuator  versus  gap 
between  the  surfaces. 
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Figure  SI.  Actuation  characteristics  of  paraffin  wax,  from  [3].  (a)  Volume  versus  temperature  of 
hexatriacontane  (98%)  under  varying  applied  stresses,  measured  by  high-pressure  dilatometry. 
(b)  Comparison  of  energy  density  for  widely  used  active  materials. 
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Figure  S2.  CNT-paraffin  film  that  was  infiltrated  with  paraffin  then  cooled  without  maintaining 
the  applied  compressive  load.  The  SEM  image  shows  large  cracks  that  prevent  vertical  actuation 
of  the  film  upon  subsequent  heating,  because  the  segments  of  the  composite  instead  expand 
laterally  while  the  paraffin  melts  and  expands. 
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Figure  S3.  Video  (100X  actual  speed)  showing  the  side  view  of  a  CNT-paraffin  nanocomposite 
actuator  upon  three  cycles  from  room  temperature  (approximately  20°C)  to  175°C. 
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Figure  S4.  Strain-temperature  curves  of  CNT-paraffin  actuators  that  were  fabricated  under 
different  amounts  of  compressive  strain  (relative  to  their  initial  height). 
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Figure  S5.  Characterization  of  self-activating  CNT-paraffin  thermal  switch  versus  temperature: 
calculated  composite  thermal  conductivity,  combined  thermal  conductivity  of  composite  and  air 
gap  thermal  conductivity,  and  measured  air  gap. 
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Figure  S6.  Performance  of  various  actuator  types,  from  Huber  et  al.  3.  (a)  Approximate  ranges 
for  the  characteristics  and  (b)  upper  limits  of  actuation  stress  versus  actuation  strain  for  various 
actuators,  (c)  Available  data  for  actuation  stress  versus  actuation  strain  for  a  CNT-paraffin  film 
actuator. 
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Appendix:  Thermal  resistance  network  model  used  to  calculate  CNT-paraffin  thermal 
conductivity 


The  thermal  conductivity  of  the  CNT-paraffin  composite  was  calculated  from  a  thermal 
resistance  network  model,  which  is  based  on  the  experimental  setup  shown  in  Figure  4A  and  is 
drawn  below. 
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In  this  model,  Qin  and  Qout  respectively  denote  the  heat  flux  into  and  out  of  the  system,  Q  latent 
is  the  heat  absorbed  by  the  composite  to  raise  its  temperature,  Thot  and  Tcold  respectively  denote 
the  temperature  of  the  heat  source  and  the  heat  sink,  RCOnvection  denotes  the  convective  thermal 
resistance  of  the  air,  Rcontact  denotes  the  contact  resistance,  Rsi  denotes  the  resistance  of  the 
silicon  substrate,  RCOmposite  denotes  the  resistance  of  the  CNT-paraffin  composite,  and  Rgap 
denotes  the  resistance  of  the  air  gap.  Based  on  the  model,  the  total  thermal  resistance  is  then 
calculated  using 
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Furthermore,  the  total  thermal  resistance  can  be  calculated  using  the  measured  heat  flux  into  the 
heat  sink  and  the  temperature  difference  between  the  heat  source  and  the  heat  sink: 


Rtotal 


Thot  TCoid 
Qout 


(2) 


Combining  Equations  1  and  2  allows  us  to  solve  for  R composite,  since  R convection  and  R contact 
have  been  independently  measured  and  the  values  for  Rsi  and  Rgap  (corresponding  to  air 
thermal  resistance)  are  known.  Last,  the  thermal  conductivity  of  the  CNT-paraffin  composite  is 
calculated  by  taking  the  inverse  of  RCOmposite •  This  model  assumes  steady  state  conditions. 
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Chapter  2:  Replica  molding  of  liquid  crystal  polymer  microstructures  for  active  surfaces 

Contents  of  this  chapter  is  in  preparation  for  publication  as:  H.  Zhao,  J.J.  Wie,  D.  Copic,  A.W. 
Orbaek,  S.  Kim  and  A.J.  Hart,  “Replica  molding  of  liquid  crystal  polymer  microstructures  for 
active  surfaces.  ” 

Introduction 

Surfaces  with  dynamic  micro-  and  nanoscale  textures  can  potentially  be  used  to  manipulate 
friction  drag,  optical  reflection,  thermal  transport,  and  other  properties.  However,  in  order  to 
realize  these  applications,  surfaces  must  be  made  consistently  over  large  areas,  and  integrated 
with  existing  laminate  and  skin  materials  for  use  in  vehicles  such  as  aircraft. 

Liquid  crystal  polymers  are  interesting  candidate  materials  for  active  surfaces  due  to  their 
tunable  mechanical  properties  1_3,  anisotropic  expansion  upon  heating  4,  and  opportunity  for 
optically  directed  actuation  5"7.  For  example,  thin  sheets  made  from  stiff  glassy  liquid  crystal 
network  (LCN)  polymers  have  been  shown  to  oscillate  at  30Hz  with  a  tip  displacement 
amplitude  of  170°  8.  However,  to  our  knowledge,  miniaturization  of  LCN  actuators  in  the  form 
of  surface-bound  micropillars  has  not  been  demonstrated. 

We  present  a  replica  molding  (RM)  fabrication  technique  9  to  obtain  high  fidelity  surface 
bound  glassy  LCN  microstructure  actuators  over  large  areas.  The  removal  of  oxygen  during 
curing  was  critical  to  achieve  high-fidelity  replicas  with  smooth  surface  texture,  while  the 
necessary  LCN  alignment  for  actuation  was  obtained  via  curing  in  a  magnetic  field.  Cast 
microstructures  had  lateral  dimensions  of  a  few  microns  and  aspect  ratios  (AR)  up  to  15: 1.  A 
gradient  based  edge  detection  algorithm  10  was  used  to  quantify  the  LCN  microstructure 
actuation,  verifying  anisotropy  of  the  network. 

Methods 

LCN  microstructure  arrays  were  fabricated  by  replica  molding  within  a  magnetic  field  in  a 
controlled  atmosphere.  The  experimental  setup  includes  a  custom  magnetic  circuit  designed  to 
straighten  and  confine  a  magnetic  through  the  sample.  The  magnetic  yoke  was  placed  inside  a 
vacuum  jar,  whose  aluminum  base  was  in  contact  with  a  hotplate  and  is  shown  in  Fig.  la.  A 
green  LED  light  source,  used  for  cross-linking,  was  located  outside  of  the  vacuum  jar. 
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A  LC  powder  mixture  composed  of  20  wt%  optically  active  4,40-bis[6-(acryoloxy)hexyloxy] 
azobenzene  monomers  (A6Z6  from  BEAM  Co.),  78  wt  %  thermally  active  RMM491  monomer 
mixture  (EMD  Millipore),  and  2  wt  %  1-784  photoinitiator  (Ciba)  was  created  The  constituents 
of  the  RMM491  mixture  and  the  optically  active  LC  monomer,  labeled  A6Z6,  are  shown  in  Fig. 
lb. 

Master  templates  can  be  fabricated  using  standard  silicon  dry  etching  processes,  or 
alternatively  CNT/SU-8  composite  microstructures  may  be  used  to  fabricate  complex  3D 
geometries  9.  From  the  master,  a  polydimethylsiloxane  (PDMS)  negative  was  cast  and  used  to 
create  the  LCN  replicas.  PDMS  was  vacuum  cast  over  the  deep  reactive  ion  etch  (DRIE) 
patterned  silicon  master  and  cured  at  120°C  for  20min  to  form  the  negative  mold  (Fig.  lc).  The 
powder  mixture  was  added  to  the  negative  and  heated  to  110°C  for  10  min  at  a  pressure  below  1 
Torr  (Fig.  Id).  The  magnetic  field  (0.3T)  applied  by  the  fixture  causes  benzene  containing  LC 
monomers  to  align  in  a  cooperative  fashion  along  it.  The  cooperative  alignment  of  benzene  rings 
inside  a  magnetic  field  was  first  suggested  by  de  Gennes  11  and  successfully  demonstrated  for 
use  in  liquid  crystal  polymers  by  Buguin  et  al.  4  The  mixture  was  cooled  to  75°C  at  l°C/min 
resulting  in  a  LC  in  the  nematic  phase  with  its  director  aligned  along  the  magnetic  field.  Placing 
the  sample  in  a  strong  magnetic  field,  which  has  been  straightened  via  a  magnetic  yoke,  and 
cooling  it  slowly  achieves  the  necessary  order  and  alignment  for  actuation.  The  sample  was  then 
exposed  to  green  light  (540  nm  wavelength)  for  40  min,  which  initiates  cross-linking  of  the 
ordered  LC  monomers  into  a  polymer  network.  After  cooling  to  room  temperature  the  LCN 
replica  was  manually  demolded  from  the  soft  PDMS  negative  (Fig  le). 

A  custom-built  microscope  was  used  to  image  the  microstructures  during  actuation  (Fig.  If). 
This  microscope  has  both  a  linearly  polarized  blue  light  and  a  white  light  source.  During 
actuation  the  samples  were  exposed  to  linearly  polarized  blue  light  (-440  nm  wavelength) 
aligned  parallel  to  the  LC  director  for  3-5  minutes.  After  actuation  the  sample  was  exposed  to 
visible  light  for  15  minutes  or  heated  to  90°C  causing  it  to  relax  and  recover  its  original  shape. 
Azobenzene  moieties  predominantly  take  on  the  trans  isomerization  during  the  relaxation  and 
the  polymer  network  recovers  its  initial  higher  order  nematic  phase. 

Results  and  Discussion 
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We  found  that  removal  of  oxygen  was  crucial  for  achieving  high  fidelity  LCN  replicas  with 
low  surface  roughness.  Samples  cast  without  degassing  have  large  voids  due  to  air  bubbles 
trapped  in  the  molten  powder  mixture  and  have  very  large  surface  roughness  (Fig.  2a).  Fig.  2b 
shows  a  sample  that  was  vacuum  cast  at  1  Torr  to  remove  air  bubbles  from  the  molten  powder 
mixture,  after  which  air  was  reintroduced  during  curing.  The  bottom  part  of  the  SEM  image 
shows  the  sample  edge  that  was  exposed  to  the  air.  Microstructures  closer  to  the  exposed  edge 
have  significantly  larger  surface  roughness,  suggesting  that  oxygen  exposure  impacts  the 
crystallization  of  the  LCN.  We  found  that  samples  with  large  surface  roughness  have 
significantly  reduced  photogenerated  strains. 

Using  these  findings  to  drive  process  improvements,  we  show  (Fig.  3)  cast  microstructures 
with  a  variety  of  sizes  and  shapes,  with  both  small  and  large  aspect  ratios.  Large  microstructures 
with  different  cross-section  shapes  are  shown  in  Lig.  3a,  while  large  arrays  of  high  aspect  ratio 
(AR,  height:width)  structures  are  shown  in  Lig.  3b.  The  semi-cylindrical  micropillars  are  55  pm 
tall  and  have  an  AR  of  approximately  15:1. 

The  LCN  micropillars  can  be  actuated  by  polarized  blue  light  capable  of  providing  sufficient 
energy  to  drive  both  the  trans-cis  and  the  cis-trans  isomerization  (Lig.  4a),  due  to  an  overlap 
between  the  cis  and  trans  isomer  absorption  spectra  12.  The  azobenzene  trails  isomer  moieties 
absorb  significantly  less  polarized  light  if  they  are  oriented  perpendicular  the  polarization, 
gradually  creating  a  statistical  buildup  of  such  isomers.  This  buildup  disturbs  the  nematic  phase, 
reducing  its  order  and  achieving  a  volumetric  shape  change  (Lig.  4b).  This  shape  change  can  be 
used  to  alter  either  the  height  of  a  micropillar  if  the  director  is  oriented  along  its  axis  (Lig.  4c)  or 
the  pillar  cross  section  if  the  director  is  oriented  perpendicular  to  it.  Because  the  actuation  is 
based  on  the  disruption  of  the  nematic  LC  phase,  having  good  alignment  prior  to  cross-linking 
the  LC  monomers  is  crucial.  A  comparison  between  similar  LCN  compositions  and  other 
actuating  materials  is  shown  in  Lig.  4d  13.  Although  the  maximum  generated  stress  of  LCN 
materials  are  lower  than  that  of  other  comparable  polymers,  we  show  that  LCN  materials  can  be 
cast  into  a  verity  of  useful  microstructures. 

An  interpolation-based  edge  tracking  algorithm  14  was  implemented  to  resolve  sub-pixel 
actuation  of  cast  LCN  micropillars,  imaged  using  the  microscope  described  above.  An  SEM 
image  and  a  top-down  optical  image  of  the  optically  active  LCN  microstructure  are  shown  in 
Lig.  5a  and  5b,  respectively.  We  measured  a  photogenerated  strain  of  0.25%  orthogonal  to  the 
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LCN  director  and  a  0.11%  contraction  along  it  after  5  min  of  exposure  to  linearly  polarized  blue 
light  oriented  perpendicular  to  the  LC  director  (Fig.  5c).  The  contraction  along  the  director 
occurs  because  the  rod-like  stiff  LC  moieties  gradually  rotate  away  from  it.  Samples  cured 
without  the  magnetic  field  showed  no  actuation.  Previous  studies  with  similar  LCN 
compositions  have  measured  a  photogenerated  strain  along  the  director  of  0.15%  and  a  stiffness 
of  1.8  GPa  using  dynamic  mechanical  analysis  (DMA)  15.  We  expect  that  larger  strains  could  be 
achieved  via  use  of  UV  light,  by  achieving  improved  LCN  alignment  during  the  replica  molding 
process,  by  composition  reformulation  at  the  expense  of  reduced  stiffness  ',  and  by  physical 
aging  of  the  material  16. 


Conclusion 

In  conclusion,  we  have  fabricated  and  tested  LCN  microactuators,  exhibiting  up  to  0.25% 
strain  along  the  director  and  a  0.11%  contraction  orthogonal  to  it  after  5  min  exposure  of 
polarized  blue  light.  For  the  first  time,  centimeter-scale  areas  of  high  aspect  ratio 
microstructures  were  fabricated,  enabled  by  careful  environmental  process  control  and  removal 
of  oxygen  during  casting.  With  further  characterization  and  improvement  of  the  optically- 
induced  strain,  this  technology  would  be  suitable  for  fabrication  of  active  surface  textures,  such 
as  for  use  in  microfluidics,  and  possibly  on  the  surfaces  of  miniature  vehicles. 
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Fig.  1.  (a)  Experimental  setup  used  to  cast  LCN  micro-structures,  (b)  Casting  of  the  PDMS 
negative  from  the  DRIE  etched  silicon  master,  (c)  Close  up  of  molten  LC  mixture  inside  the 
magnetic  field,  (d)  Manual  demolding  of  replicated  structures,  (e)  Diagrams  of  the  monomers 
used  in  the  LCN  mixture,  (f)  Schematic  of  custom  microscope  schematic  used  during  testing. 
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Figure  2.  (a)  LCN  replicas  cast  without  degassing  of  the  molten  LC  mixture,  (b)  Degassed 


replica  that  was  exposed  to  air  during  curing  with  close  up  images  of  selected  regions,  showing 
increased  surface  roughness  in  the  presence  of  oxygen. 
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Figure  3.  (a)  Heterogeneous  low  aspect  ratio  optically  active  LCN  microstructures,  (b)  High 
yield,  large  area  array  of  asymmetric  micro  structures.  Inset  shows  close  up  image  of 
microstructures. 
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Figure  4.  (a)  trcins-cis-trans  isomerization  of  azobenzene  due  to  polarized  blue  light  exposure, 
(b)  Schematic  of  nematic  and  isotropic  LCN  phase  containing  thermally  active  (purple)  and 
optically  active  LC  molecules  (orange),  (c)  Volumetric  shape  change  due  to  LCN  phase  change, 
(d)  Maximum  stress  and  strain  generated  by  comparable  actuating  materials  13. 
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Figure  5.  (a)  SEM  image  of  optically  active  LCN  micropillar,  (b)  Sample  optical  image  of  the 
top-down  view  from  a  custom  microscope  used  in  for  edge  tracking.  Arrows  indicate  strain 
direction,  (c)  Results  of  gradient  based  edge  tracking  algorithm,  showing  edge  movement. 
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Chapter  3:  Strain-engineered  manufacturing  of  freeform  carbon  nanotube  microstructures 

Contents  of  this  chapter  have  been  published  as:  M.De  Voider,  S.  Park,  S.  Tawfick  and  A.J.  Hart, 
“Strain-engineered  manufacturing  of  freeform  carbon  nanotube  microstructures.  ”  Nat. 

Commun.  5,  4512  (2014). 

Introduction 

Scalable  fabrication  of  microstructures  that  mimic  the  hierarchical  surface  designs  found  in 
nature  has  been  a  long-standing  aspiration  of  material  scientists  '"5.  While  symbiotic  growth  of 
the  integrated  circuit  (IC)  and  micro  electro  mechanical  systems  (MEMS)  industries  has  enabled 
innovations  in  3D  fabrication  that  leverage  semiconductor  processing  tools,  these  methods,  such 
as  interference  or  inclined  exposure  lithography  are  typically  limited  to  arrays  of  identical 
structures  6"8.  Rapid  prototyping  methods  such  as  direct  laser  writing,  multiphoton  lithography, 
and  focused  ion  beam  milling  can  create  arbitrary  forms  but  are  serial,  and  therefore  have  lower 
areal  throughput  910.  It  is  also  especially  difficult  to  fabricate  surface  structures  having  curved 
and/or  re-entrant  geometries. 

On  the  other  hand,  use  of  locally  directed  actions,  such  as  mechanical  stresses,  capillary 
forces,  and  electromagnetic  fields,  along  with  their  interactions  with  templates,  offers 
opportunity  to  create  novel  self-organized  geometries  and  to  design  fabrication  processes  that 
achieve  attractive  combinations  of  dimensional  control  and  throughput  2-8,11'15.  Examples  abound 
in  soft  materials  and  chemical  systems  including  micro-scale  reaction-diffusion  patterns  16,17, 
self-assembly  of  block  copolymers  18,  and  helical  aggregation  of  polymer  nanopillars  19-20. 
However,  many  of  these  processes  need  further  development  to  achieve  structural  uniformity 
over  large  substrates. 

We  present  a  novel  approach  taking  advantage  of  microscale  top-down  lithographic 
patterning  in  conjunction  with  nanoscale  self-organization,  which  enables  large-area  fabrication 
of  freeform  micro  structures  made  of  aligned  carbon  nanotubes  (CNTs).  This  process  leverages 
the  influence  of  the  catalyst-substrate  interactions  on  the  growth  rate  of  CNTs  21 ,  creating  stress 
gradients  during  synthesis  that  guide  the  CNTs  into  curved  microscale  geometries.  Our  method 
is  analogous  to  the  well-known  use  of  thin  film  stress  to  create  curved  and  folded  MEMS 
structures  22.  However,  because  our  process  is  based  on  an  additive  chemical  synthesis  process 
instead  of  a  subtractive  etching  and  release  technique,  it  enables  the  direct  synthesis  of  complex 
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microstructures  that  are  perpendicular  rather  than  parallel  to  the  substrate.  This  has  two  major 
implications:  it  enables  fabrication  of  closely  packed  arrays  of  structures  with  heterogeneous 
shapes,  and  the  porosity  of  the  CNT  forests  enables  conformal  coating  after  growth  to  modify 
chemical  and/or  mechanical  properties.  We  demonstrate  this  latter  point  by  conformal  coating  of 
CNT  “microtruss”  arrays  by  atomic  layer  deposition  (ALD)  and  polymer  CVD,  which  increases 
their  mechanical  stiffness  without  changing  the  geometry. 

Methods 

Substrate  patterning  and  CNT  growth:  The  catalyst  and  TiN  layers  are  patterned  on  (100)  silicon 
wafers  with  300  nm  of  thermally  grown  SiOi.  Each  layer  is  patterned  by  lift-off  processing,  by 
photolithography  (photoresist  1X845)  followed  by  ultrasonic  agitation  in  acetone.  The  TiN  layer 
is  deposited  and  patterned  first,  and  then  the  catalyst  layer  (1  nm  Fe  upon  10  nm  AI2O3)  is 
deposited  and  patterned.  The  wafer  is  then  diced  into  ~lxl  cm  pieces,  and  the  substrates  are 
placed  in  the  quartz  tube  furnace,  and  the  CNT  growth  is  performed.  The  recipe  starts  by  with 
flowing  100/400  seem  of  He/tb  while  heating  to  775  °C  over  10  minutes  (ramping  step);  then 
the  system  is  held  at  775  °C  for  10  minutes  (annealing  step)  while  maintaining  the  gas  flow. 

Then  100  seem  of  C2H4  is  added  to  the  gas  mixture  at  775  °C  for  CNT  growth  for  the  desired 
duration.  The  typical  growth  rate  is  ~50  pm/minute  on  Fe/AbCE/SiCF.  Once  the  CNTs  have 
grown,  C2H4  is  removed  from  the  gas  mixture  and  the  furnace  is  cooled  to  <100  °C.  After 
cooling,  the  system  is  purged  with  He  before  the  sample  is  removed.  Optionally,  C2H4  flow  can 
be  maintained  while  cooling  down  to  improve  the  adhesion  of  the  CNT  microstructures  to  the 
substrate.  Once  the  cooling  step  is  complete,  the  quartz  tube  is  purged  with  1000  seem  of  He  for 
5  minutes  before  opening  up  the  end  caps  and  retrieving  the  samples. 

Small  Angle  X-ray  Scattering  (SAXS):  SAXS  characterization  of  CNT  microstructures  was 
performed  at  the  Cornell  High  Energy  Synchrotron  Source  (CHESS)  using  the  G1  beamline 
(10+0.1  keV,  0.13nm  wavelength).  The  beam  is  focused  to  a  10pm  spot  using  a  single  bounce 
monocapillary.  The  CNT  sample  is  placed  on  a  motorized  stage  and  the  focused  X-ray  beam  is 
passed  through  the  sample.  The  scattered  beam  is  collected  using  a  2D  detector  and  the  measured 
intensities  were  normalized  to  the  original  intensity  measured  by  another  detector  at  upstream  of 
the  CNT  sample.  The  scattering  data  is  then  fitted  to  a  mathematical  model  assuming  a  log- 
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normal  distribution  of  hollow  cylinders  to  calculate  the  CNT  diameters  as  well  as  the  Herman’s 
parameter  for  CNT  alignment.  The  detailed  procedure  is  described  by  Bedewy  et  al.  30. 

Results  and  Discussion 

Fabrication.  CNTs  grown  by  chemical  vapour  deposition  (CVD)  from  a  high-density 
arrangement  of  catalyst  nanoparticles  on  a  substrate  are  known  to  self-organize  into  vertically 
aligned  assemblies  often  called  “CNT  forests”  23,24.  We  first  observed  that  the  density  and  rate 
of  CNT  forest  growth  from  a  widely  studied  supported  catalyst  (Fe/AbCB,  1/10  nm)  can  be 
influenced  by  the  material  immediately  beneath  the  catalyst.  This  premise  is  shown  in  Fig.  la; 
patterning  of  CNT  growth  catalyst  (Fc/AbCb)  on  a  SiCb/TiN  “checkerboard”  followed  by 
exposure  to  standard  CVD  conditions  (see  Methods)  results  in  a  “bi-level”  CNT  micropillar 
array.  The  catalyst  patterns  directly  on  SiCb  grow  CNTs  to  -100  pm  (in  <  2  minutes),  whereas 
the  patterns  on  TiN  (upon  SiCb)  grow  CNTs  to  50  pm  in  the  same  time  span.  In  Fig.  lb,  “tri¬ 
level”  CNT  forests  are  grown  by  arranging  patches  of  catalyst  on  SiCb,  70  nm  TiN,  and  140  nm 
TiN.  This  principle  could  be  extended  to  an  arbitrary  number  of  levels  or  even  continuous  height 
gradients  via  additional  lithography  and  underlayer  deposition  steps  that  modulate  the  growth 
rate  via  catalyst-substrate  interactions. 

Next,  we  used  the  above  differential  growth  principle  to  design  compound 
catalyst/underlayer  patterns  that  directly  form  curved  CNT  forest  geometries.  If  a  continuous 
micro-scale  catalyst  pattern  is  placed  partially  on  SiCb  and  partially  on  TiN,  the  differential 
growth  rates  induce  stress  within  the  CNT  microstructure.  For  example,  as  shown  in  Fig.  2a  and 
2b,  a  square  catalyst  pattern  with  half  of  its  area  on  the  TiN  layer  bends  toward  the  side  which  is 
upon  TiN,  due  to  the  difference  in  growth  rate  on  the  coupled  halves  of  the  structure.  The  stress 
is  transferred  between  contacting  CNTs  at  the  boundary  region  via  mechanical  entanglement  and 
van  der  Waals  interactions  among  the  CNTs.  Depending  on  the  curvature  and  length  of  the 
structures,  slanted  micropillars  (Fig.  2a),  or  arches  (Fig.  2b)  can  be  fabricated.  Because  the  local 
interaction  and  differential  growth  rate  determines  the  trajectory  of  each  structure,  large  arrays 
with  nearly  identical  anisotropic  shapes  can  be  produced  as  shown  in  the  SEM  images. 
Importantly,  these  3D  structures  are  fabricated  using  only  two  standard  photolithography  steps, 
one  for  patterning  the  TiN  layer,  and  one  for  patterning  the  catalyst  layer. 
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The  curvature  can  be  controlled  by  designing  the  amount  of  overlap  between  the  catalyst  and 
the  TiN  underlayer.  This  is  illustrated  in  Fig.  3a  and  b,  which  respectively  show  arrays  of  round 
and  square  cross-section  micropillars  where  the  overlap  distance  is  varied  from  left  to  right 
(increments  of  5  pm).  As  expected,  the  portion  of  the  pillars  growing  on  TiN  is  always  shorter, 
and  as  a  result,  all  pillars  bend  towards  the  TiN  side.  As  the  portion  of  overlap  decreases,  the 
stress  induced  by  the  differential  growth  rate  causes  increased  bending  (smaller  radius  of 
curvature),  reaching  a  maximum  when  the  catalyst  shape  is  split  symmetrically  by  the  TiN  layer. 
With  <50%  overlap  on  TiN,  the  curvature  increases  gradually  until  the  structure  is  only  slightly 
curved  at  the  rightmost  extent  of  the  array.  The  CNTs  are  generally  tangential  to  the  curvature  of 
the  micro  structures,  similar  to  the  CNT  alignment  observed  in  CNT  forests. 


Static  model  of  stress-driven  CNT  curvature.  The  coupling  of  stress  and  CNT  growth  rate, 
via  the  anisotropic  mechanics  of  the  CNT  forest  25,  is  a  complex  problem.  However,  we  find  that 
a  static  mechanical  model  analogous  to  that  used  for  differential  expansion  of  a  bi-material 
cantilever  beam  can  reasonably  predict  the  curvature  of  CNT  structures  grown  from  an 
overlapped  catalyst/TiN  rectangle.  Starting  from  the  classical  formulation  of  the  bimetallic  strip 
model  26,  we  replace  the  temperature-dependent  expansion  term  by  a  differential  lengthening 
term  representing  the  CNT  growth  rate.  Accordingly,  the  curvature  of  the  compound  CNT 
microstructure  is  described  as: 


i 

P  w^3(l+m)2+(l+mn)(mz  +^))  V  ' 

Here  p  is  the  radius  of  curvature,  R1  and  Rz  are  the  growth  rates  (1  denotes  CNTs  on 
Fe/AECE/SiOi  and  2  denotes  CNTs  on  Fe//Al2C>3/TiN),  and  w  is  the  CNT  micropillar  width.  In 


addition,  m  and  n  are  defined  as 
m  =  —  ,n  —  —  (2) 

Wj  E1  v  ' 

where  w  denotes  the  width  and  E  denotes  the  respective  Young's  Moduli  of  the  segments  .  The 


value  of  n  is  specified  as  0.6  which  is  the  ratio  of  the  measured  areal  mass  density  of  CNTs  on 
the  respective  underlayers;  however,  because  the  elongation  of  each  layer  is  specified  in  the 
model,  the  output  is  insensitive  to  this  value.  The  geometric  parameters  are  defined  in  Fig.  3c. 
Using  the  calculated  curvatures,  and  the  weighted  average  growth  rate,  the  shapes  of  the 


resultant  CNT  microstructures  were  visualized  using  Matlab.  The  simulation  results  correspond 
to  the  rows  of  structures  in  the  SEM  image  Fig.  3b.  To  compare  the  experiments  to  the 
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simulation,  the  tip  position  was  characterized  in  horizontal  and  vertical  axes,  normalized  to  the 
base  dimension  (w),  as  shown  in  Fig.  3c.  For  both  the  experiment  and  simulation,  the  x  position 
of  the  tip  reaches  its  maximum  at  0.4  overlap,  and  the  y  position  reaches  its  minimum  at 
approximately  0. 6-0.7  overlap. 

The  differences  between  the  predicted  and  calculated  displacements  arise  because  the  model 
does  not  capture  the  exact  kinetics  of  CNT  growth,  which  varies  with  time.  Moreover,  it  cannot 
consider  how  the  stress  between  the  two  portions  of  the  structure,  which  are  idealized  as 
perfectly  coupled  without  slip,  influences  the  deformation.  CNT  forests  have  anisotropic 
mechanical  properties,  with  the  lateral  stiffness  (perpendicular  to  the  CNT  alignment)  typically 
much  less  than  the  axial  stiffness  27 ,  therefore  in  principle  favoring  greater  deflection  due  to 
built-in  stress  gradient.  Local  wrinkling  and  buckling  of  the  CNTs  in  the  compound 
microstructures  indicates  that  the  growth  stress  causes  complex  mechanical  deformations,  which 
cannot  be  predicted  by  linear  elasticity  and  the  bimaterial  deflection  model  (see  Fig.  4a). 

To  gain  further  insight  into  the  mechanical  coupling  causing  stress-driven  bending,  we 
designed  a  “striped”  structure  where  alternate  catalyst/underlayer  regions  are  coupled  with  a 
large  interfacial  area.  This  structure  is  symmetric,  so  it  grows  straight  vertically  yet  has 
significant  internal  stresses.  As  shown  in  Fig.  4a  and  the  neighboring  insets,  the  faster  growing 
CNTs  deform  collectively  into  a  wavelike  pattern.  Therefore,  while  the  vertical  growth  rate  of 
the  structure  is  matched  at  the  interface,  the  faster  growing  side  still  accumulates  longer  CNTs 
and  these  CNTs  bend  and  possibly  buckle  to  accommodate  their  additional  length.  This 
deformation  mode  is  similar  to  what  is  observed  in  mechanically  compressed  forests  28. 

As  seen  in  Fig.  3a  and  supplementary  figure  la,  some  structures  separate  at  the  interface 
between  the  differentially  growing  regions,  due  to  shear  stresses  at  the  interface.  In  future  work, 
this  could  be  remedied  by  adapting  the  structure  design  to  reduce  the  interfacial  shear  stress 
between  the  segments,  or  by  using  a  more  gradual  growth  rate  gradient  to  accommodate  the 
stress  gradient,  such  as  a  three-layer  design  as  in  Fig.  lb.  However,  we  importantly  find  that 
structures  that  do  not  separate  during  synthesis  can  withstand  large  subsequent  deformations 
without  failure.  For  instance,  supplementary  figure  lb  shows  an  SEM  image  of  an  arrangement 
of  bent  pillars  during  compression  to  50%  vertical  strain,  at  which  point  delamination  at  the 
interface  begins  only  where  the  pillars  kink  near  their  midsection. 
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CNT  and  catalyst  morphology.  The  strain-engineered  CNT  microstructures  show  differing 
CNT  density  and  alignment  in  the  fast  and  slow  growing  portions.  In  Fig.  4a,b,  the  CNTs  grown 
from  catalyst  on  TiN  appear  to  have  greater  vertical  alignment  influenced  by  the  interface  with 
the  faster- growing  region.  On  the  other  hand,  the  CNTs  grown  from  catalyst  on  Si02  are  less 
aligned,  due  to  the  retarding  force  from  the  slower-growing  mating  regions. 

We  hypothesized  that  the  differential  CNT  growth  behavior  on  the  TiN  underlayer  may  be 
attributed  to  differences  in  the  catalyst  morphology,  which  can  influence  the  CNT  diameter 
distribution  and  number  density.  Atomic  Force  Microscopy  (AFM)  showed  that  as-deposited 
catalyst  layers  on  SiOi  and  TiN  have  a  similar  topology  (supplementary  figure  2).  However, 
upon  annealing  in  H2/He  prior  to  hydrocarbon  exposure  (see  methods),  the  TiN-catalyst  layer 
forms  shallow  mounds,  tens  of  nanometers  high  and  hundreds  of  nanometers  wide,  in  addition  to 
smaller  catalyst  particles.  The  control  case  of  catalyst  on  SiCF  does  not  exhibit  such  topography. 
For  samples  with  80  nm  TiN  layer,  the  average  catalyst  particle  height  and  spacing  were 
calculated  to  be  5.2  nm  and  19  nm  respectively.  Compared  to  those  on  SiCF  (7.5  nm  and  18  nm, 
respectively  29),  the  catalyst  particle  sizes  are  smaller  on  average  while  the  spacing  is 
comparable.  In  addition,  the  Root  Mean  Square  (RMS)  roughness  of  the  annealed  catalyst/TiN 
layer  is  5.6  nm,  which  is  considerably  higher  than  that  of  catalyst  on  SiCF  (1.4  nm). 

Small  Angle  X-ray  Scattering  (SAXS)  was  used  to  further  investigate  the  CNT  forest 
morphology  30,3 '.  For  this  experiment,  CNTs  were  grown  for  10  minutes  on  SiCF,  and  on  40  nm, 
and  80  nm  TiN  layers;  these  samples  reached  lengths  of  800,  500,  and  400  pm  respectively.  The 
scattered  X-ray  intensities  were  fitted  to  a  mathematical  form  factor  model  for  hollow  cylinders 
32  to  calculate  the  diameter  and  Herman’s  orientation  parameter,  which  is  a  measure  of 
alignment.  Both  the  CNT  diameter  (Fig.  4e)  and  alignment  (supplementary  figure  3)  are  shown 
to  be  lesser  for  increasing  TiN  underlayer  thickness.  Specifically,  CNTs  on  SiCF  have  initial 
average  diameter  of  9.5  nm,  while  CNTs  on  TiN  are  approximately  8  nm  in  diameter.  These 
measurements  further  support  the  AFM  data  which  showed  that  the  catalyst  particles  on  TiN 
layers  are  smaller  than  those  on  SiCF.  The  decrease  in  diameter  with  growth  time  has  been 
attributed  to  diffusion  of  the  catalyst  into  the  AI2O3  30,33 .  The  Herman’s  orientation  parameter 
increases  from  the  top  of  the  forest  (the  initial  growth),  then  reaches  a  maximum,  and  then 
decreases  toward  the  bottom  of  the  forest.  This  trend  is  typically  observed  for  CNT  forests 
grown  by  CVD,  and  has  been  attributed  to  density  variation  due  to  collective  activation  and 
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deactivation  of  the  growing  CNT  population  30.  The  measured  areal  mass  density  of  the  CNT 
forests  is  0.01 1  mg  mm"2  on  80  nm  TiN  and  0.018  mg  mm"2  on  SiCF.  Therefore,  these  methods 
consistency  show  that  placement  of  the  TiN  layer  under  the  catalyst  results  in  CNT  forests  with  a 
smaller  average  CNT  diameter,  lesser  alignment  of  CNTs,  and  a  lower  density.  Although  we 
have  not  directly  compared  the  mechanical  properties  of  the  different  segments,  we  expect  that 
the  CNT  forests  on  TiN  have  lower  stiffness,  and  lesser  anisotropy  in  mechanical  properties. 

The  positive  correlation  between  CNT  density  and  alignment  is  also  expected  based  on  recent 
literature,  which  demonstrated  that  CNT-CNT  crowding  controls  alignment  within  non-patterned 
forests  34 . 

Based  on  our  understanding  of  the  elementary  catalyst/underlayer  designs  that  achieve 
unidirectional  bending,  we  designed  a  variety  of  more  complex  patterns  that  produce  exemplary 
CNT  microstructures  having  complex  curvature.  For  instance,  a  compound  shape  consisting  of  a 
“+”  catalyst  microfeature  with  each  arm  offset  by  a  rectangular  TiN  underlayer  results  in  growth 
of  twisted  CNT  microstructures  (Fig.  5a),  resembling  macroscale  propellers.  The  first-order 
applicability  of  the  bimaterial  bending  model  discussed  above  inspired  us  to  evaluate  the 
suitability  of  finite  element  modeling  (FEM)  to  predict  the  shapes  of  these  structures.  These 
were  simulated  using  Comsol  FEM  software  as  illustrated  in  Fig.  5a,  capturing  the  uncoupled 
differential  growth  rate  (Fig.  la)  as  a  50%  expansion  mismatch,  and  estimating  the  Young’s 
modulus  as  30  MPa  for  the  CNTs  on  TiN  and  50  MPa  on  SiCF  35. 

Similarly,  thin  semicircles  of  CNTs  can  be  directed  to  curve  outward  by  offsetting  the  TiN 
underlayer  as  shown  in  Fig.  5b.  Further  structural  complexity  is  shown  by  the  scrolling  of  thin 
offset  rectangular  patterns  (Fig.  5c).  Last,  exotic  hierarchical  arrangements  can  be  formed  by  the 
interaction  of  closely  spaced  structures,  such  as  the  self-organization  of  offset  circular 
micropillars  into  wavy  patterns  (Fig.  5d)  that  are  reminiscent  of  macroscale  crochet  stitching. 

We  hypothesize  that,  after  the  individual  structures  bend  unidirectionally  and  contact  one 
another,  their  continued  growth  and  steric  hindrance  causes  the  wavy  pattern  to  form.  More 
investigation  is  needed  to  understand  the  complex  deformations  of  these  structures,  and  their 
relationship  to  the  mechanics  of  the  CNT  forest  and  the  mechanical  feedback  on  the  growth 
process  itself. 

Notably,  in  spite  of  the  complex  geometries  and  local  deformations,  all  of  these  structures 
can  be  produced  with  impressive  consistency  over  large  arrays.  Arrays  of  several  hundred 
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structures  were  examined  and  shown  to  exhibit  nearly  identical  forms,  with  defects  most 
frequently  arising  from  debris  due  to  the  lithography  process  rather  than  the  CNT  growth  step.  In 
this  study,  we  explored  structures  with  critical  dimensions  as  small  as  5  pm  (Fig.  5),  and  found 
that  the  uniformity  of  the  structures  was  not  sacrificed  at  this  scale.  We  expect  that  smaller  3D 
microstructures  could  be  made  while  still  using  optical  lithography  along  with  high-precision 
alignment  of  the  catalyst  and  TiN  layers.  Notably,  sub-micron  vertical  CNT  features  have  been 
fabricated  for  use  as  interconnects  36.  Further,  because  the  curved  and  twisted  geometries  result 
from  collective  behavior  of  CNT  forest  growth,  we  expect  the  structures  to  require  a  certain 
minimal  size  to  average  out  the  stochastic  variation  of  individual  CNT  growth  rate  and  catalyst 
lifetime.  Considering  the  approximate  CNT-CNT  spacing  within  the  current  microstructures 
(-100  nm),  the  minimal  feature  size  of  reliable  growth  of  curved  microstructures  may  be  limited 
to  ~1  pm,  though  this  requires  further  investigation. 

Last,  we  show  that  3D  CNT  structures  can  be  post-processed  via  both  wet  and  dry  methods 
that  enable  tuning  of  their  properties  and  functionality.  Low-density  bent  CNT  micropillars  can 
be  transformed  into  robust  densely  packed  CNT  structures  by  capillary  forming  (Fig.  6a-b) 

35.37.38  y0  do  So5  the  substrate  is  exposed  to  a  stream  of  heated  acetone  vapor,  causing  acetone  to 
condense  onto  the  CNTs  and  substrate,  and  infiltrate  each  CNT  microstructure.  Upon  subsequent 
evaporation  of  the  acetone,  the  CNT  forest  shrinks  laterally,  due  to  the  surface  tension  of  the 
shrinking  meniscus.  Previously,  we  showed  that  capillary  forming  of  vertical  CNT 
microstructures  increases  the  Young’s  modulus  in  compression  approximately  100-fold,  from 
~50  MPa  to  4  GPa  35.  These  values  are  comparable  to  soft  rubbers  and  stiff  epoxies, 
respectively.  Moreover  in  the  present  case  the  capillary  forming  process  preserves  the  curved 
geometry  and  increases  the  lateral  deflection. 

Alternatively,  curved  CNT  microstructures  can  be  coated  conformally  via  vapor  phase 
methods,  thereby  enabling  decoupled  control  of  geometry  and  mechanical  properties.  To 
investigate  this,  we  fabricated  CNT  “microtruss  arrays”  (Fig.  6c-e),  which  are  analogous  to  truss 
designs  used  in  composite  materials  to  achieve  high  strength  and  energy  absorption  at  relatively 
low  density  39-41 .  The  CNT  microtrusses  each  consist  of  four  corner  members  and  a  central 
pillar,  meeting  at  an  apex.  We  explored  coating  of  the  CNT  microtrusses  with  both  parylene  (by 
chemical  vapor  deposition,  CVD,  supplementary  figure  4)  and  alumina  (by  atomic  layer 
deposition,  ALD,  Fig.  6d).  Upon  vapor  phase  infiltration  of  the  precursors,  the  CNTs  and 
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bundles  within  the  forest  are  coated  individually  and  conformally  42 ,  enabling  fine-tuning  of  their 
porosity  and  mechanics  without  altering  the  microstructure  geometry.  Via  flat  punch 
compression  testing,  we  found  that  a  51  nm  AI2O3  coating  on  the  CNTs  increases  the  mechanical 
stiffness  by  more  than  100-fold;  typical  loading-unloading  curves  are  shown  in  Fig.  6e  and 
supplementary  figure  4  for  different  coating  materials  and  thicknesses.  The  equivalent  stiffness 
range  of  the  3D  CNT  microtrusses  is  0.36  to  54  kN  m-1,  which  spans  typical  values  of  MEMS 
springs  used  in  probe  card  arrays  43.  A  further  attribute  of  the  TiN  underlayer  is  its  electrical 
conductivity.  In  the  future,  electrical  integration  of  strain-engineered  CNT  structures  in 
conjunction  with  the  post-processing  methods  described  above  could  be  useful  in  advanced 
microsystems,  including  as  structural  elements  or  microsensors  44,45.  Such  applications  could 
take  advantage  of  the  thermal  and  mechanical  durability  of  CNTs,  as  well  as  the  anisotropic 
properties  arising  from  their  alignment  and  collective  curvature. 

The  ability  to  fabricate  large  arrays  of  3D  microstructures  is  also  conducive  to  mimicry  of 
nature’s  advanced  functional  materials  5.  We  say  this  because  the  structures  shown  in  Fig.  2  and 
5  have  similar  size  and  geometry  to  butterfly  wings  that  have  anisotropic  wetting  properties  46, 
dry  adhesive  contacts  on  the  legs  of  beetles  and  gecko  lizards  3,  and  microscale  sensing  hairs 
found  on  spiders  47.  Therefore,  we  suggest  that  3D  CNT  microstructures  offer  opportunities  to 
further  engineer  the  excellent  dry  adhesive  48  49  and  superhydrophobic  50  properties  previously 
reported  for  CNT  forests.  For  mechanical  applications  such  as  dry  adhesion,  an  important 
consideration  is  the  adhesion  of  the  structures  to  the  substrate,  as  well  as  the  mechanical 
resilience  of  the  structures  upon  repeated  loading.  While  we  have  found  that  the  segments  of  the 
curved  structures  are  well  adhered  to  each  other  (supplementary  figure  lb),  we  anticipate  that 
coating  or  transfer  steps  may  be  necessary  to  anchor  the  structures  for  use  as  robust  surface 
contacts  49 .  Exploration  of  potential  optical  and  photonic  applications  would  require  further 
miniaturization  of  the  structures.  Nevertheless,  combined  with  the  emerging  methods  to  grow 
and  pattern  CNT  forests  on  large-area  substrates  51"53,  we  are  optimistic  that  the  methods  shown 
here  can  enable  large-area  surface  coatings  having  advanced  functionalities. 

Conclusion 

In  summary,  we  show  that  strain-engineered  CNT  growth  enables  the  scalable  fabrication  of 
complex  3D  microstructures  having  unidirectional  and  multidirectional  curvature.  The 
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capability  to  produce  such  structures  en  masse  using  only  2D  patterning  methods  along  with 
standard  thermal  processing  contrasts  the  limitations  of  many  existing  processes  that  require 
serial  processing  or  sequential  exposure  using  complex  inclined  lithography  methods.  Moreover, 
the  structures  can  profit  from  the  mechanical  robustness  and  electrical  conductivity  of  CNTs,  and 
their  mechanical  and  surface  properties  can  be  engineered  independently  from  their  geometry  by 
conformal  coating  of  the  CNTs  after  growth.  This  represents  a  highly  attractive  principle  for 
materials  design,  and  is  promising  for  scalable  manufacturing  of  3D  microstructured  surfaces 
having  biomimetic  properties. 
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Figure  1.  Multi-height  CNT  micropillar  growth  using  TiN  as  a  rate-controlling  catalyst 
underlayer,  (a)  Process  steps  including  two  lithography  steps,  and  accompanying  (color  added) 
pattern  of  cylindrical  CNT  micropillars  grown  on  “checkerboard”  of  alternating  TiN  (80nm 
thickness)  squares  on  SiCF  on  the  right  in  a  black  box.  Scale  bar  is  40  pm  (b)  Triple-height  arrays 
made  by  substrate  and  catalyst  patterning  with  three  lithography  (TiN  layers  are  40  nm  and  80 
nm  thickness).  Scale  bars  are  500pm  (left)  and  200  pm  (right). 
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Figure  2.  Curved  CNT  micropillars  grown  from  catalyst  rectangles  partially  overlapping  TiN 
(shown  as  green  layer  in  the  schematics).  The  pillars  bend  toward  the  TiN  side,  which  grows 
more  slowly  and  couples  to  the  faster-growing  region  of  catalyst  directly  on  SiCF.  (a)  Closely- 
spaced  short  structures  with  -45°  takeoff  angle.  Scale  bars  are  100  pm  (panel)  and  20  pm  (inset), 
(b)  Arch-like  structures  that  curve  over  and  contact  the  substrate  at  their  distal  ends.  .  Scale  bars 
are  200  pm  (panel)  and  20  pm  (inset). 
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Figure  3.  Shape  gradients  and  their  use  in  elementary  modeling  of  the  stress-driven  CNT  growth 
process.  (a,b)  arrays  of  round  and  square  CNT  micropillars  with  decreasing  catalyst/TiN  overlap 
from  left  to  right.  Scale  bars  are  100  pm.  (c)  Schematic  of  geometric  parameters  used  to  quantify 
deformation  of  square  micropillars,  (d)  Output  of  differential  expansion  model  predicting 
vertical  and  lateral  deflection  versus  substrate  pattern  design,  comparing  simulations,  inset  to  (b), 
and  measurements  from  SEM  images. 
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Figure  4.  Analysis  of  CNT  morphology,  (a)  SEM  image  of  micropillar  with  catalyst  on  fine 
stripes  of  TiN,  with  close-up  images  showing  the  CNT  morphology  difference  and  evidence  of 
internal  stress  causing  local  deformation.  Scale  bars  are  100  pm  (left),  2  pm  (upper  right)  and  10 
pm  (lower  right),  (b)  Close-up  images  of  the  aligned  CNT  morphology  in  adjacent  segments  of 
the  micropillar.  Scale  bars  is  500  nm.  (c,d)  SEM  images  of  the  substrate  morphology  after 
annealing  in  H2  and  rapid  cooling,  showing  hierarchical  morphology  of  nanoparticles  on  the 
catalyst/TiN  area.  Scale  bars  are  500  nm.  (e)  Diameter  mapping  within  the  CNT  forests,  from 
fitting  of  SAXS  data.  Red  -  CNT  of  Fe/AbC>3,  Green  -  CNT  on  Fe/AhCE/TiN  (40nm),  Blue  - 
CNT  on  Fe/AbOs/TiN  (80nm). 
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Figure  5.  Design  and  fabrication  of  arrays  of  complex  microstructures.  Each  of  the  four  panels 
contains  a  sketch  (top  left)  of  the  catalyst-TiN  pattern  used  for  CNT  growth,  a  close-up  SEM 
image  (top  right),  and  and  an  exemplary  array  of  structures  (bottom),  (a)  Twisted  propeller- like 
structures,  made  from  azimuthally  offset  crosshair  catalyst/TiN  layers,  also  compared  to  finite 
element  model  (FEM)  prediction  of  shape  (top  right).  Scale  bars  are  40  pm  (top  right),  and  500 
pm  (bottom),  (b)  Outward  curving  semicircles  made  from  radially  offset  catalyst/TiN  layers. 
Scale  bars  are  5  pm  (top  right),  and  100  pm  (bottom),  (c)  Scroll-like  deformation  of  thin- walled 
microstructures  with  slight  catalyst/TiN  overlap.  Scale  bars  are  10  pm  (top  right),  and  250  pm 
(bottom),  (d)  Collective  organization  of  bending  microstructures  into  a  wavy  pattern.  Scale  bars 
are  10  pm  (top  right),  and  40  pm  (bottom). 
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Figure  6.  Post-processing  of  3DCNT  microstructures.  (a,b)  Array  of  bending  micropillars  before 
and  after  capillary  forming,  respectively.  The  unidirectional  anisotropic  morphology  is 
maintained  and  the  micropillar  cross  section  decreases,  while  the  lateral  deflection  angle  of  the 
structures  increases.  Scale  bars  are  200  pm.  (c)  Microtruss  array  fabricated  to  study  mechanical 
reinforcement  by  vapor-phase  coating  of  the  CNTs.  Scale  bars  are  100  pm  (left)  and  400  pm 
(right),  (d)  SEM  images  of  aligned  CNTs  in  the  forest  sidewall  before  and  after  coating  with 
AI2O3  by  ALD.  Scale  bars  are  1  pm.  (e)  Load-displacement  curves  before  and  after  ALD 
coating,  where  the  stiffness  is  measured  from  the  slope  of  the  unloading  curve.  Green  -  No 
coating,  Blue  -  AI2O3  coating  (ALD,  100  cycles),  Red  -  AI2O3  coating  (ALD,  1000  cycles). 
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Figure  S 1 .  SEM  image  of  a  CNT  microstructure  where  the  fast  (on  Fe/AbCb/SiCb)  and  slow  (on 
Fe/AbCb/TiN/SiCb)  growing  portions  delaminated  completely  at  the  vertical  interface,  (b)  In-situ 
SEM  indentation  of  a  CNT  micro-truss.  The  curved  beams  do  not  split  apart  even  after  a  large 
deformation,  except  locally  near  the  center  where  the  large  deformation  causes  kinking. 


Figure  S2.  AFM  images  comparing  topography  of  annealed  TiN  layer,  and  Fe/AbCb  catalyst  on 
TiN.  Images  show  that  TiN  layer  roughens  when  annealed  with  features  that  are  100s  of 
nanometers  wide  and  10s  of  nanometers  high  (red).  Annealed  catalyst  over  TiN  layer  shows 
catalyst  particles  on  larger  scale  mounds  (blue).  The  insets  show  AFM  scans  before  annealing. 
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Figure  S3.  Herman's  parameter  mapping  within  the  CNT  forests,  from  fitting  of  SAXS  data. 
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Figure  S4.  Parylene-coated  CNT  microtrusses  for  mechanical  reinforcement,  (a)  Overview  and 
close  up  SEM  images  of  the  coated  CNT  trusses,  (b)  load-displacement  curve  of  a  single 
microtruss  with  identical  geometry  as  Ah03-coated  trusses  shown  in  the  main  text. 
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Chapter  4:  Scaling  the  stiffness,  strength,  and  toughness  of  ceramic-coated  nanotube  foams 
into  the  structural  regime 

Contents  of  this  chapter  have  been  published  as:  A.  Brieland-Shoultz,  S.  Tawfick,  S.J.  Park,  M. 
Bedewy,  M.  Maschmann,  J.W.  Baur  andA.J.  Hart,  “Scaling  the  stiffness,  strength,  and 
toughness  of  ceramic-coated  nanotube  foams  into  the  structural  regime.  "Adv.  Funct.  Mater.  24, 
5728-5735  (2014). 

Introduction 

Design  of  the  hierarchical  structure  and  composition  of  cellular  materials  offers  a  versatile 
means  to  tune  mechanical  properties,  such  as  Young’s  modulus,  strength,  and  toughness. 

Cellular  materials  can  also  achieve  unique  combinations  of  properties,  while  maintaining  large 
surface  areas  and  low  mass  density.  For  example,  the  mechanical  properties  of  metallic  foams 
enable  their  use  for  energy  absorption  in  consumer  goods  packaging  and  composite  sandwich 
structures,  and  their  porosity  enables  use  as  battery  and  capacitor  electrode  supports  '.  Ceramic 
foams  are  widely  used  as  mechanical,  thermal  and  acoustic  insulators,  for  flow  filtration  and 
regulation,  and  as  catalyst  supports  2.  Many  natural  materials  use  analogous  design  principles; 
for  example,  the  exquisite  structure-function  relationship  of  bone  spans  from  stiff  and  dense 
cortical  load  bearing  bone  to  highly  porous  cancellous  bone  where  haematopoiesis  occurs  3. 
Further,  the  geometric  and  compositional  hierarchy  of  soft  and  hard  constituents  result  in  the 
high  toughness  of  bone  at  low  mass  density  4. 

Recently,  micro-  and  nanofabrication  methods  such  as  photo,  X-ray,  and  interference 
lithography  have  been  applied  to  prototype  micro-  and  nanoscale  foams  with  unique  cell 
geometries  and  hierarchical  structures  5.  For  example,  fabrication  of  epoxy  lattices  by 
interference  lithography  enabled  a  priori  tuning  of  the  cell  size  in  the  sub-micron  regime.  These 
foams  showed  a  classical  power-law  relationship  between  Young’s  modulus  and  relative  density 
(E  ~  \p/ps\n),  with  7?=1.26;  at  pr  =  0.5  the  energy  absorption  reached  4  J/cm3  at  strain  <  0.25  5. 

Use  of  sacrificial  polymer  templates  to  make  hollow  tube  metallic  microlattices  resulted  in  a 
foam  with  record  low  density  (pr  =  p/ps  =  2x10  4,  p  ~1  mg/cm3)  and  high  recoverability,  but  low 
Young’s  modulus  (E  ~  550  kPa)  6.  In  this  case,  due  to  the  hollow  structure  of  the  metal  truss 
members,  the  relationship  between  modulus  and  relative  density  followed  a  power  law  with  n= 2. 
This  contrasts  most  fibrous  nanoscale  aerogels  and  foams,  having  77=2-3,  which  is  attributed  to 
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low  connectivity  between  the  constituent  struts  1 .  When  comparing  foam  materials,  a  higher 
power  law  exponent  results  in  a  relatively  larger  change  in  E  for  the  same  change  in  density. 

While  these  methods  have  enabled  exploration  of  new  design  and  property  regimes  for 
engineered  foams,  their  manufacturing  cost  and  complexity  currently  prohibit  them  from  scaling 
to  large  quantities  that  can  be  incorporated  in  macroscale  engineering  systems.  In  contrast,  most 
industrially  manufactured  foams  have  isotropic  and  random  cell  geometry.  Closed  cell  metallic 
and  polymer  foams  are  usually  manufactured  by  gas  injection  into  melts,  while  open  cell  metallic 
foams  are  usually  made  by  vapor  deposition  or  electrodeposition  of  metal  precursors  on  polymer 
scaffolds.  Porous  metallic  and  ceramic  materials  can  also  be  manufactured  by  compaction  and 
sintering  of  powders  8. 

Directed  growth  and  self-organization  of  nanostructures  is  emerging  as  an  attractive  and 
potentially  scalable  route  to  create  hierarchical  materials  having  both  composition-  and 
geometry-derived  mechanical  properties.  Carbon  nanotubes  (CNTs)  have  outstanding  properties 
including  high  Young’s  modulus,  strength,  conductivity,  and  thermal  stability.  Moreover,  CNTs 
can  be  manufactured  in  large  quantities  by  cost-effective  chemical  vapor  deposition  (CVD) 
methods  9,  and  assembled,  either  during  or  after  CVD,  into  ordered  forms  such  as  vertically 
aligned  “forests”,  yarns,  and  sheets  10.  While  as-grown  CNT  forests  have  high  compressibility  n- 
12  and  have  shown  promise  as  energy-dissipating  surfaces  and  interlayers  u’ 13,  CNT  forests  and 
isotropic  foams  made  from  CNTs  13  have  low  density  usually  in  the  range  of  10-100  kg/m3  or 
p/ps  =  0.5-5%  (considering  ps  to  be  equal  to  solid  graphite  at  2200  kg/m3). 

Therefore,  it  remains  an  opportunity  to  explore  how  self-assembly  of  CNTs,  in  combination 
with  conformal  coating  methods  14,  can  be  used  to  create  bulk  materials  with  efficient  and  widely 
tunable  properties,  and  to  compare  these  properties  to  commercially  available  structural 
materials.  In  this  study,  we  show  that  nanoscale  coating  of  a  ceramic  onto  clean  CNT  forests 
enables  unprecedented  tuning  of  Young’s  modulus,  strength,  and  toughness,  without  changing 
the  structure  or  composition  of  the  nanoporous  material. 

Methods 

CNT  growth:  The  catalyst  for  CNT  forest  growth  is  deposited  on  a  (100)  silicon  wafer  with  300 
nm  of  thermally  grown  silicon  dioxide.  Photoresist  (SPR220)  is  patterned  by  photolithography 
and  developed,  and  then  films  of  10  nm  AI2O3  and  1  nm  Fe  are  deposited  by  sputtering  (Kurt  J. 
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Lesker  PVD  Lab  18).  The  photoresist  is  removed  by  performing  two  consecutive  cycles  (8 
minutes  each)  of  ultrasonication  in  acetone,  followed  by  flushing  with  IPA.  CNT  growth  is 
performed  by  thermal  CVD  in  a  quartz  tube  furnace  (Thermo-Fisher  Minimite,  22  mm  inner 
diameter).  Pieces  of  the  patterned  catalyst  wafer  are  placed  in  the  quartz  tube  and  sealed,  and 
then  the  furnace  is  flushed  with  1000  seem  He  for  5  minutes.  Then  100/400sccm  H2/He  is  flowed 
for  20  minutes;  the  temperature  is  ramped  up  to  775°C  during  the  first  10  minutes  and  held  for 
the  next  10  minutes.  For  CNT  growth,  100  seem  C2H4  is  added  to  the  gas  mixture  while  the 
temperature  is  maintained  at  775°C.  Once  the  growth  step  is  complete,  the  furnace  power  is 
turned  off,  and  the  tube  is  flushed  with  1000  seem  He  until  the  furnace  is  below  100°C.  The 
height  of  the  CNT  microstructures  was  controlled  by  the  growth  time,  with  the  growth  rate 
equaling  approximately  50  pm/minute. 

ALD  coating:  After  growth,  CNT  micropillars  were  coated  with  alumina  deposited  via  atomic 
layer  deposition  (ALD).  The  ALD  process  is  performed  in  a  low  pressure  CVD  chamber  (Oxford 
OpAl)  at  300  mTorr  and  150°C  via  alternating  pulses  of  A1(CH3)3  and  H2O  with  a  purge  between 
the  pulses. 

Mechanical  testing:  Micropillars  and  perforated  forests  were  tested  in  compression  using  a  MTS 
Nanoindenter  XP  with  flat  diamond  punch  (100  pm  diameter).  In  situ  SEM  compression  tests 
were  performed  using  a  flat  nanoindentation  tip  (100  pm  diameter)  within  a  custom-built 
micromechanical  test  frame  inside  a  SEM  (FEI  QUANTA  ESEM),  as  described  by  Maschmann 
et  al  1 1 . 

X-ray  scattering:  As-grown  CNT  micropillars  as  well  as  ALCh-coatcd  micropillars  were 
mounted  on  a  motorized  sample  stage  in  the  Synchrotron  X-ray  beamline.  Energy  of  the  X-ray 
beam  used  was  10  +/-  0.1  keV.  Scattering  patterns  were  collected  on  a  2D  area  detector.  The 
size  of  the  X-ray  beam  was  controlled  using  a  pair  of  motorized  precision  slits.  More  details  of 
the  experimental  setup  and  the  data  analysis  techniques  are  published  in  Bedewy  et  al  19. 

Results  and  Discussion 

Fabrication  and  compression  testing  of  ceramic-carbon  nanotube  foams 

First,  CNT  forests,  comprising  vertically  aligned  multi-wall  CNTs,  were  synthesized  from  a 
supported  metal  catalyst  film  (see  Methods)  by  atmospheric  pressure  chemical  vapor  deposition 
(CVD)  15.  Then,  AI2O3  was  deposited  by  atomic  layer  deposition  (ALD),  conformally  coating 
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the  CNTs  and  CNT  bundles  within  the  forest,  resulting  in  a  nanoporous  composite  with  ceramic- 
carbon  core-shell  structural  members.  The  thickness  of  AI2O3  on  the  CNTs  was  controlled  by 
the  number  of  ALD  process  cycles.  During  each  cycle  of  ALD,  the  precursors  (A1(CH3)3  and 
H2O)  are  sequentially  introduced  in  two  half-reaction  stages,  and  the  thickness  during  each  cycle 
is  self-limiting  to  approx.  1.3  A  . 

The  change  in  morphology  of  the  CNTs  due  to  AI2O3  coating  is  shown  in  Fig.  la-c.  First, 
ALD  AI2O3  nucleation  occurs  preferentially  at  surface  defects  on  the  CNTs  (Fig.  SI1)  17  and  at 
contact  points  among  CNTs.  After  nucleation,  the  coating  grows  to  fully  cover  the  individual 
CNTs,  and  in  some  cases  bundles  of  CNTs.  After  a  certain  thickness  the  coated  CNTs  and 
bundles  coalesce.  Judging  from  SEM  images  of  the  forest  sidewalls,  this  coalescence  begins  at 
coating  thickness  of  approximately  30  nm,  compared  to  the  average  CNT-CNT  spacing  of 
approximately  50  nm  (number  density  ~3xl010/cm2).  The  local  CNT  spacing  and  density 
represent  a  distribution,  determined  by  the  waviness  and  bundling  of  the  CNTs  which  can  be 
approximated  by  their  average  orientation  angle  as  discussed  later. 

CNT  micropillars,  made  using  a  lithographically  patterned  catalyst,  were  used  to  study  the 
scaling  of  mechanical  properties  with  the  number  of  ALD  cycles,  and  to  relate  the  mechanics  to 
the  penetration  depth  of  ALD  coating.  Test  samples  were  fabricated  with  arrays  of  CNT 
micropillars  having  radii  ranging  from  10-35  pm  (Pig.  Id).  The  micropillars  were  tested  in 
quasi-static  axial  compression  using  a  nanoindenter  with  a  cylindrical  flat  diamond  tip  (50  pm 
radius).  The  Young’s  modulus  was  obtained  from  the  unloading  slope  in  continuous  modulus 
mode  (10  nm  oscillations  at  50  Hz)  at  a  compressive  pre-strain  of  0.05  (approximately  2  pm 
compression  of  the  micropillar). 

Load-displacement  curves  of  individual  15pm  radius  AI2O3-CNT  micropillars  with  30  pm 
height,  with  coatings  of  0,  100,  and  1000  ALD  cycles,  are  compared  in  Lig.  le.  The  load 
capacity  of  the  pillars  increases  from  <0.1  mN  for  the  uncoated  (bare  CNT)  pillar,  to  0.7  mN 
with  13  nm  coating  (100  ALD  cycles),  and  130  mN  with  51  nm  coating  (1000  cycles).  As 
shown  in  Lig.  2c  for  Rpiiiar  =  15  pm,  the  Young’s  modulus  scales  from  14.5  MPa  for  uncoated 
CNTs  to  24.5  GPa  with  51  nm  coating. 

Structural  characterization  of  coating  morphology 
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To  further  understand  the  structure-property  relationship  via  ALD  coating,  we  measured  the 
coating  thickness  and  the  density  evolution  per  ALD  cycle.  To  determine  the  penetration  depth, 
we  coated  patterned  CNT  forest  “line”  features  and  found,  as  shown  in  Fig.  3a,  that  the 
penetration  depth  is  approximately  5-7  pm.  To  overcome  this  limitation  and  demonstrate  the 
scaling  on  large  sample  areas,  perforated  “honeycomb”  patterns  (Fig.  3b)  with  10  pm  spacing 
between  holes  were  synthesized.  The  perforations  enabled  uniform  ALD  coating  throughout  the 
CNT  forest,  and  therefore  were  used  for  subsequent  characterization. 

Transmission  SAXS  was  then  used  to  characterize  the  morphology  of  the  as  grown  and 
AbCb-coated  CNTs  by  placing  the  honeycomb  film  sample  in  the  beamline  path  of  a  synchrotron 
X-ray  beam  and  collecting  the  scattered  X-rays  on  an  area  detector  (Fig.  4,  S2)  18- 19.  As  in  our 
previous  work  19 ,  the  statistical  distribution  of  outer  diameters  for  both  coated  and  uncoated 
CNTs  was  obtained  by  fitting  a  mathematical  model  of  polydisperse  core-shell  cylinders  to  line- 
scans  of  scattering  intensity.  This  determines  the  average  outer  and  inner  radii  of  the  bare  CNTs 
as  5.0  and  3.5  nm,  respectively.  The  linescans  are  obtained  by  integration  of  the  scattered 
intensities  within  a  sector  of  +/-  10°  about  the  axis  of  maximum  intensity  on  the  2D  scattering 
pattern,  as  shown  in  Fig.  S2.  As  seen  in  Fig.  4a,  the  peak  corresponding  to  the  scattering  from 
the  population  of  AbCb-coated  CNTs  shifts  to  a  lower  value  on  the  x-axis  with  increasing 
coating  thickness.  This  indicates  a  significant  increase  in  the  size  of  the  scatters,  which  is 
quantified  by  fitting  the  data  for  coated  samples  with  a  core-shell  scattering  model.  Analysis  of 
the  coated  samples  gives  average  coating  thickness  of  8.5,  14.5  and  26.7  nm  for  50,  100  and  200 
cycles  of  ALD,  (Fig.  4b).  For  100  and  200  cycles,  a  second  peak  on  the  SAXS  linescans  is  seen 
at  a  q  value  higher  than  the  form-factor  peak.  This  peak  could  arise  from  the  observed 
roughening  of  the  AI2O3  layer  at  the  greater  coating  thickness  (Fig.  S2c). 

We  also  estimated  the  average  coating  thickness  by  measuring  the  weight  of  the  samples 
before  and  after  ALD  coating.  We  found  a  sub-linear  relationship  between  mass  and  the  number 
of  cycles,  showing  that  the  deposited  thickness  per  ALD  cycle  is  affected  by  the  ability  of  the 
precursors  to  diffuse  into  the  CNT  forest;  this  reduces  as  the  porosity  of  the  sample  decreases 
when  subject  to  multiple  cycles.  We  therefore  constructed  a  model  where  the  pore  size  is 
decreased  with  each  subsequent  cycle.  We  found  (Fig.  4b)  that  the  best  fit  to  the  SAXS  and 
weight  measurements  corresponds  to  a  decrease  of  effective  diffusion  coefficient  of  vapor  phase 
precursors  in  porous  media  by  0.5%  per  cycle,  starting  at  the  140th  cycle  of  ALD  coating.  We 
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did  not  consider  the  increase  in  available  area  due  to  the  increase  in  effective  diameter,  which  is 
counteracted  by  the  coalescence  of  coated  CNTs  as  the  coating  thickness  increases.  Improved 
uniformity  of  coating  could  likely  be  achieved  by  increasing  the  precursor  cycle  duration. 

As  a  result  of  diffusion-limited  coating,  the  Young’s  modulus  of  the  AI2O3-CNT  composites 
depended  strongly  on  the  micropillar  diameter.  The  full  matrix  of  compression  tests  performed 
with  pillars  having  10-35  pm  radius,  at  each  coating  thickness,  is  summarized  in  Fig.  S3.  Here, 
modulus  is  inversely  proportional  to  pillar  diameter,  as  the  reinforcement  of  larger  pillars  is 
limited  by  ALD  penetration  and  therefore  the  final  composite  pillar  is  a  reinforced  shell 
surrounding  a  soft  core  of  uncoated  CNTs.  These  trends,  and  specifically  the  size-dependent 
modulus  of  non-coated  pillars,  are  also  attributed  to  size-dependent  growth  effects  on  the  density 
of  CNTs  within  the  forest.  This  topic  is  currently  under  separate  study  in  our  group  20. 

Model  of  cellular  mechanics 

We  now  develop  a  model  of  the  Young’s  modulus  of  the  AI2O3-CNT  composites  under 
quasi-static  loading.  First,  as  verified  by  SAXS,  we  note  that  ALD  coating  creates  a  nanoporous 
composite  where  the  load-bearing  elements  are  nanoscale  core-shell  beams.  These  beams  are 
loaded  primarily  in  bending  when  the  structure  is  subject  to  compressive  loading.  Previous  work 
using  in  situ  SEM  digital  image  correlation  has  shown  that  the  CNT-CNT  contact  points  within 
an  uncoated  CNT  forest  do  not  slip  upon  compressive  loading  to  strains  below  the  failure 
threshold  of  approximately  5%  local  strain  21 .  Asa  result,  uncoated  CNT  forests  can  be 
approximated  as  open-cell  foams  where  the  deformation  of  individual  CNTs  and  bundles 
represent  the  cell  walls. 

We  approximate  the  nanoporous  network  of  interconnected  wavy  AI2O3-CNT  beams  as  an 
effective  cellular  solid  with  rectangular  unit  cell  geometry  (Fig.  2a).  We  consider  the  CNTs  as 
an  array  of  beams  loaded  simultaneously  in  bending  and  compression  22,  having  identical 
diameter  and  length.  We  further  use  SAXS  to  determine  the  average  orientation  angle  of  each 
CNT  with  respect  to  the  vertical  axis  (9  =  40°,  with  a  random  distribution  of  azimuthal  angles)  23, 
giving  the  unit  cell  aspect  ratio  (HcNT/WCeii=  1.67).  Although  the  CNT  alignment  and  therefore 
the  orientation  are  known  to  vary  throughout  the  forest  due  to  changes  in  the  CNT  density,  here 
we  use  a  single  average  value  to  represent  the  effective  properties. 
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With  this  information,  we  first  calculate  the  effective  unit  cell  size,  where  each  unit  cell 


contains  one  angled  CNT  having  length  (Lcnt)  and  inclination  angle  (0).  This  is  an 
approximation  for  the  real  structure  of  the  forest,  which  certainly  contains  a  distribution  of  CNT 
angles,  effective  lengths,  and  effective  unit  cell  sizes.  The  stiffness  of  a  unit  cell  kcNT  is 
represented  by  a  cantilever  beam  in  combined  axial  compression  and  end-loaded  cantilever  beam 
bending  (Fig.  2a, b), 

1  _  sin2  6  3  cos2  6 

-  -  —  7  Lcnt  +  —  7  Lcnt 

(Eq.  1) 
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Here,  Ecnt  is  the  Young’s  modulus  (1  TPa  24,  based  on  the  wall  area  only),  and  the  cross- 
sectional  area  of  the  CNT  ( Acnt )  and  second  moment  of  area  (/cat)  are  both  determined  from  the 
average  CNT  outer  and  inner  diameters  measured  by  SAXS.  The  total  stiffness  of  the  foam  is 
therefore  related  to  the  unit  cell  stiffness  by  the  assembly  of  M  (across  the  horizontal  cross- 
section)  by  N  (along  the  height)  unit  cells  (Fig.  2a)  into  a  cylindrical  solid  representing  a 
micropillar, 
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We  derive  Lcnt  =  210  nm,  as  well  as  N  and  M,  which  give  h  =  160  nm  and  w  =  95  nm. 

We  now  predict  the  modulus  of  the  composite  assuming  uniform  coating  of  the  slanted 
CNTs  within  each  unit  cell,  causing  a  proportional  contribution  of  each  cell  to  the  overall 
stiffness  of  the  micropillar  structure.  Therefore,  we  then  scale  the  Young’s  modulus  of  the  unit 
cell,  using  a  core-shell  beam  model  where  the  AI2O3  shell  thickness  is  determined  by  the  number 
of  ALD  cycles.  The  AI2O3  Young’s  modulus  is  taken  to  be  180  GPa  25. 


(El), 


—  f1  T  +  J7  T 

composite  ^ CNT 1  CNT  ^  ^ AI2Q31  A/203  * 


(Eq.  3) 


We  used  ECNT  equal  to  The  comparison  between  the  measured  compressive  modulus  and 
the  predicted  modulus  from  the  unit  cell  model  is  shown  in  Fig.  2c.  For  coating  thickness  below 
~5  nm,  the  effective  spring  model  with  constant  Lcnt  closely  predicts  the  experimentally 
measured  scaling  of  elastic  modulus  with  coating  thickness.  As  a  result,  we  conclude  that  the 
primary  mechanism  of  axial  deformation  of  the  foam  is  bending  of  the  individual  nanoscale 
elements.  However,  to  our  surprise,  assuming  constant  Lcnt  does  not  accurately  predict  the 
measured  stiffness  at  coating  thickness  >10  nm.  This  is  because  the  length  of  the  slanted  core- 
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shell  beams  decreases  as  more  coating  is  deposited,  causing  additional  junctions  to  form  between 
nearby  CNTs  whose  gaps  are  filled  by  the  thicker  coating.  To  adapt  the  model  to  the 
morphology  evolution,  we  iteratively  update  the  characteristic  length, 

LcNT\n+1  —  LCNT\„  —  2  tn+1 

Here,  n  and  n+1  are  the  cycle  numbers,  and  1  is  the  AI2O3  deposition  thickness  during  the 
cycle,  which  is  determined  from  the  SAXS  and  density  analysis.  This  is  schematically  shown  in 
Fig.  2d.  The  variable  Lcnt  model  accurately  predicts  the  modulus  scaling  at  large  coating 
thickness,  up  to  the  maximum  coating  of  5 1  nm.  At  this  point,  the  increased  reinforcement  due  to 
coalescence  of  individual  CNTs  and  bundles  (as  seen  in  Fig.  lc)  changes  the  effective  unit  cell 
height,  thus  increasing  load  transfer  by  compression  of  the  structural  members  rather  than 
bending.  Moreover,  further  increase  in  density  of  the  composite  is  limited,  as  we  suspect  the 
exterior  voids  are  fully  closed  and  the  coating  precursors  can  no  longer  penetrate  uniformly  and 
deeply  into  the  foam.  As  a  result,  we  calculated  the  Young’s  modulus  to  represent  ALD 
penetration  depth  of  5  pm  from  the  outside  surfaces  of  the  micropillar.  In  Fig.  S4  we  provide 
plots  showing  the  dependence  of  Young’s  modulus  on  the  ALD  penetration  depth. 

While  penetration  of  larger  micro-scale  geometries  is  limited  under  the  present  ALD 
conditions,  the  perforated  honeycomb  design  (Fig.  3c)  shows  potential  to  achieve  uniform 
composite  properties  over  large  areas,  as  it  enables  uniform  infiltration  of  the  coating  precursors. 
A  honeycomb  CNT  forest  with  1000  cycles  ALD  AI2O3  has  E  =  38.3  ±  3.7  GPa  based  on  the 
solid  area  of  the  cross-section.  Nevertheless,  future  study  should  focus  on  optimizing  the  depth 
and  uniformity  of  infiltration  and  coating,  as  determined  by  the  deposition  temperature,  pressure, 
and  the  precursor  chemistry. 

The  AI2O3-CNT  foams  also  withstand  large  deformations  at  high  stresses  before  failure,  and 
therefore  have  exceptional  toughness.  In  Fig.  5  we  show  the  measured  stress  versus  energy 
absorbed  for  foams  with  selected  coating  thicknesses.  The  data  here  is  obtained  by  integrating 
the  area  under  the  quasi-static  stress-strain  curves  of  micropillars  compressed  to  large  strains. 

We  note  three  regions,  as  compressive  strain  increases:  (1)  a  linear  regime  where  the  energy 
absorbed  is  proportional  to  the  applied  stress;  (2)  a  steep  increase  in  energy  absorbed, 
corresponding  to  the  plateau  in  the  stress-strain  curve;  and  (3)  densification  as  the  material  fails. 
Large  energy  absorption,  characterized  by  an  increase  in  total  energy  stored  due  to  deformation 
without  much  increase  in  the  resultant  stress,  occurs  until  the  vertical  tip  of  each  of  these  curves 
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before  the  slope  is  decreased  due  to  the  small  resultant  strain  at  the  foam  densification  regime 
(3).  This  indicates  that  most  of  the  energy  absorbed  is  due  to  the  large  deformations  of  the 
individual  cells,  a  favorable  property  for  structural  foams. 

The  high  toughness  indicates  that  the  CNTs,  even  at  very  low  weight  %  (~2.5  wt%  CNTs 
after  1000  cycles  ALD),  bridge  cracks  in  the  AI2O3  foam  that  form  upon  compressive  loading. 
This  prevents  catastrophic  fracture  until  a  very  large  strain,  and  hence  results  in  large  energy 
absorption.  We  believe  that  rapid  crack  propagation  is  mitigated  due  to  the  collective  buckling 
of  the  AI2O3-CNT  struts  held  in  proximity  by  van  der  Waals  interactions,  or  in  some  areas  by 
thin  ceramic  junctions.  In  contrast,  classical  open  cell  (e.g.,  Al)  foams  have  continuous  uniform 
struts  connected  at  the  junctions  and  spanning  the  cells.  Thus,  at  the  same  density,  classical 
foams  require  larger  forces  to  buckle  the  members,  and  exhibit  a  narrower  plateau  in  the  stress- 
strain  curve  than  observed  for  the  AI2O3-CNT  foams. 

The  deformation  and  failure  mechanisms  are  further  understood  by  examining  in  situ  SEM 
videos  of  micropillars  under  compression  (Fig.  5b).  In  these  videos  and  the  accompanying 
stress-strain  curves,  we  find  that  foams  with  approximately  10  nm  AI2O3  coating  fully  recover 
large  compressive  strains,  and  exhibit  significant  hysteresis  (i.e.,  material  damping).  We  suspect 
this  is  due  to  a  balance  between  the  reinforcement  provided  by  the  AI2O3  coating,  enabling  the 
elastic  energy  to  reversibly  overcome  the  stiction  of  the  struts  due  to  van  der  Waals  forces,  while 
remaining  flexible  to  prevent  fracture  at  large  global  strains.  At  thinner  coatings,  the  applied 
stress  is  not  recovered  as  stiction  dominates,  whereas  at  thicker  coatings,  fracture  occurs, 
resulting  in  permanent  deformation  albeit  with  very  high  toughness.  The  role  of  the  coating 
thickness  in  joining  and  bundling  CNTs  into  a  hierarchical  nanocomposite  structure  is  evident  by 
comparing  the  fracture  behavior  of  pillars  in  the  thick  coating  regime.  The  pillar  with  28  nm 
coating  splits  outward  as  fracture  propagates  parallel  to  the  CNTs,  while  the  pillar  with  43  nm 
coating  initially  cracks  then  crushes  uniformly  under  uniaxial  compression. 

In  summary  of  the  above  data,  we  achieve  tuning  of  the  elastic  modulus  by  more  than  1000- 
fold  by  simply  varying  the  thickness  of  the  AI2O3  ALD  coating.  The  achieved  values  (Fig.  5a) 
range  from  (on  the  low  end)  those  typical  of  uncoated  CNT  foams  and  other  low-density 
nano  structured  foams,  to  (on  the  high  end)  those  known  for  structural  aluminum  honeycombs 
having  millimeter- sc  ale  cell  size  26.  Moreover,  the  strength  of  our  AI2O3-CNT  composites, 
which  is  determined  as  the  maximum  stress  before  the  slope  of  the  loading  curve  decreases 
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abruptly  (Fig.  5b,  Fig.  S5),  is  10-fold  higher  than  isotropic  aluminum  foams  having  the  same 
density  of  p  =  1000  kg/m3,  and  at  this  point  is  comparable  to  compact  bone  which  has  2-fold 
higher  density  (p  =  2000  kg/m3)  3.  The  modulus  that  we  achieve  also  exceeds  that  measured  for 
composites  of  mechanically  densified  CNT  forests  with  aerospace  grade  epoxy,  where  E  =  9  GPa 
is  achieved  at  20%  CNT  volume  fraction  27 .  As  stated  before,  E  is  defined  as  the  slope  of  the 
unloading  curve,  and  crv  is  defined  as  the  location  on  the  stress-strain  curve  where  the  slope  starts 
to  decrease  after  the  initial  linear-elastic  regime. 

Comparing  the  scaling  behavior  to  conventional  open  cell  foams,  we  note  that  the  Young’s 
modulus  and  strength  of  the  AI2O3-CNT  composite  scale  as  E~ p2-8  and  Gy~p2  9  (Fig.  5a,b).  Ideal 
bend-dominated  open-cell  foams  such  as  open-cell  polymer  foams  28,  and  aluminum 
honeycombs  29  scale  as  E~p2  and  ay  ~p2.  The  modulus  scaling  of  our  ALD-coated  CNTs 
incidentally  resembles  the  Young’s  moduli  of  isotropic  CNT  foams  30,  CNT  aerogels 31,  and  Si 
aerogels  32  which  scale  as  ~p\  and  commonly  attributed  to  the  interdependence  of  strut  thickness 
and  cell  size.  This  explains  the  discrepancy  between  the  measured  E  and  the  constant  length  strut 
analytical  model  (Fig.  2c).  On  the  other  hand,  ultralight  open-cell  foams  such  as  Ni-P 
microlattices  6  show  E~p2  because  of  the  strut  aspect  ratio  is  independent  of  density.  High 
modulus  scaling  power  (n),  as  shown  by  the  AI2O3-CNT  foams,  is  advantageous  because  we 
begin  with  low-density  as  grown  CNTs  and  achieve  a  comparatively  large  increase  in 
mechanical  properties  with  increasing  density.  The  controllability  and  interdependence  of 
stiffness,  strength,  and  toughness  will  certainly  depend  on  the  initial  CNT  density,  diameter,  and 
morphology,  and  moreover  on  the  crystallinity  and  composition  of  the  coating.  Further  study  is 
warranted  to  understand  these  phenomena,  and  to  engineer  the  mechanical  energy  storage  and 
dissipation  properties 

In  the  present  study,  the  scale  of  the  materials  tested  was  limited  to  the  honeycomb  films 
covering  ~1  cm2  area,  and  this  area  is  determined  by  the  size  of  the  CVD  furnace  used  for  CNT 
growth.  Nevertheless,  the  design  principles  of  conformal  coating  and  hierarchical  patterning 
could  be  realized  on  much  larger  scales,  such  as  within  coatings  of  CNTs  grown  on  woven 
micro-fibers  or  other  porous  frameworks  33 .  Moreover,  emerging  methods  to  pattern  CNT 
growth  catalyst  using  continuous  printing  methods  such  as  laser  printing  34,  and  machines  for 
roll-to-roll  CVD  on  metal  foils,  suggest  that  large-area  coating  and  lamination  of  CNT  forests 
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will  be  possible  in  the  future.  Several  schemes  for  high-speed  continuous  ALD  have  also  been 
reported  and  are  being  commercialized  35 . 

Conclusion 

Conformal  coating  of  CNT  forests  may  serve  as  a  scalable  platform  for  the  engineering  and 
fabrication  of  hierarchical  foams  with  nanoscale  cell  geometry.  As  demonstrated  here,  AI2O3- 
CNT  foams  made  using  lithographically  patterned  CNT  microstructures  followed  by  ALD 
coating  show  1000-fold  tunable  Young’s  modulus,  compressive  strength,  and  toughness;  and  the 
scaling  of  Young’s  modulus  over  this  wide  range  can  be  accurately  predicted  by  the  mechanics 
of  beam  deformation  as  informed  by  X-ray  characterization  of  the  coating  thickness.  Moreover, 
the  AI2O3-CNT  foams  have  ~  10-fold  higher  strength  than  commercial  aluminum  foams  at  the 
same  density,  as  well  as  higher  toughness  due  to  the  high  ductility  exhibited  by  the 
nano  structured  ceramic  coating.  As  a  result,  these  foams  may  be  suitable  for  filtration  and  shock 
absorption  under  high  temperature  and  pressure  conditions.  We  believe  that  capabilities  to  grow 
and  coat  CNTs  in  hierarchically  microstructured  geometries,  on  advanced  fibers,  and  within 
porous  materials,  will  enable  realization  of  core-shell  CNT-based  foams  in  a  wide  range  of 
macroscopic  configurations.  Future  exploration  of  both  ceramic  and  metallic  coatings  could 
enable  new  high-performance  foams  for  structural  applications  as  well  as  electrochemical  energy 
storage  materials  where  porosity  and  mechanical  robustness  are  essential  to  high  performance. 
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Figure  1.  Morphology,  structure  and  compression  testing  of  A1203-carbon  nanotube  (CNT) 
foam,  (a)  Schematic  and  scanning  electron  microscope  (SEM)  image  of  as-grown  vertically 
aligned  CNTs  (“forest”).  CNT  forest  conformally  coated  with  (b)  13  mu  thick  and  (c)  51  mu 
thick  A1203  by  atomic  layer  deposition  (ALD).  (d)  Row  of  A1203/CNT  micropillars  fabricated 
for  micro-compression  testing,  (e)  Load-displacement  measurement  of  micropillars  (D  =  30  pm), 
having  nanoscale  morphology  as  shown  in  (a-c). 
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Figure  2.  Unit  cell  model  to  predict  compressive  stiffness  of  AI2O3-CNT  foams,  based  on 
bending  deformation  of  angled  core-shell  beams,  (a)  Schematic  of  composite  micropillar,  and 
representative  decomposition  into  NxM  unit  cells,  (b)  Schematic  of  an  individual  unit  cell 
showing  a  CNT  beam  in  bending,  with  load  decomposed  into  parallel  and  perpendiculat 
components,  (c)  Comparison  of  measurements  and  model  of  the  scaling  of  Young’s  modulus 
with  the  AI2O3  coating  thickness,  (d)  Schematics  of  the  evolution  of  coating  morphology, 
indicating  reduction  of  the  characteristic  beam  length  which  modifies  the  effective  unit  cell 
geometry. 
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Figure  3.  Penetration  depth  and  hierarchical  design  of  AI2O3-CNT  composite  films,  (a)  Cross 
section  SEM  of  cleaved  microwall  (100  pm  width)  showing  the  coating  penetration  depth  (bright 
areas)  from  the  top  and  the  sides  of  the  micropillar,  (b)  Cross  section  SEM  of  cleaved  microwall 
(10  pm  width)  showing  unifonn  coating  throughout  the  thickness,  (c)  Hierarchical  honeycomb 
pattern  with  10  pm  wall  thickness  designed  to  achieve  unifonn  scaling  of  ALD  coating  over 
large  areas. 
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Figure  4.  Characterization  of  core-shell  CNT-AI2O3  morphology  by  Small  Angle  X-ray 
Scattering  (SAXS)  spectra  of  “honeycomb”  pattern  films,  (a)  I-q  curves  obtained  from 
integration  of  SAXS  images,  showing  evolution  of  primary  scattering  peak  corresponding  to 
population  of  cylindrical  nanostructures,  (b)  AI2O3  coating  thickness  obtained  by  fitting  SAXS 
data  with  mathematical  model,  compared  to  measurements  derived  from  weight  change  of  each 
sample  after  ALD.  Curve  fit  to  (b)  represents  diffusion-limited  infiltration  of  the  CNT  forest  by 
the  ALD  coating  precursors. 
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1  Linear-elastic  2  Plateau  3  Densification 


Figure  5.  Toughness  of  AI2O3-CNT  foams,  (a)  Strain  energy  absorbed  versus  stress  for 
micropillars  with  different  AI2O3  coating  thickness,  showing  the  transition  from  linear  elastic 
compression  to  the  plateau  (increase  in  strain  without  increase  in  stress)  and  the  densification 
regimes  (increase  in  stress  without  increase  in  strain),  (b)  SEM  images  from  in  situ  compression 
testing  showing  the  failure  modes  for  coated  CNT  pillars  at  high  strains:  buckling  for  13  nm 
thick  coatings  and  cracking  for  28  and  42  nm  thick  coatings.  Top  row  shows  micropillars  during 
compression,  and  bottom  row  shows  the  same  respective  test  after  unloading. 
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♦  CNT-AI203  [this  study] 

CNT  foam  [Worsley  et  al.  2009] 

CNT  foam  [Worsley  et  al.  2009] 

•  Al  foam  (Alcan)  [Andrews  et  al.  1999] 


O  Al  foam  (Fraunhofer) 
[Andrews  et  al.  1999] 
AI2O3  hydrogel 
Si  aerogels 
APDMS  XSU-8 


0  Polymer-CNT  composites  [Cebeci  et  al.  2009] 
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Figure  6.  Measured  power-law  relationship  between  bulk  density  and  compressive  (a)  Young’s 
modulus  and  (b)  strength  of  the  AI2O3-CNT  foams,  compared  to  benchmark  published  and 
commercially  available  materials. 
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Figure  S 1 .  S  AXS  characterization  of  AI2O3  coating  thickness  evolution,  (a)  Schematic  of 
experimental  setup  for  high  energy  transmission  SAXS  of  vertically  aligned  CNT;  (b)  Line  scans 
of  scattering  spectra  having  various  number  of  ALD  cycles,  where  arrow  indicates  the  decrease 
in  q  with  the  increased  number  of  cycles  in  addition  to  the  appearance  of  a  new  peak  at  high  q 
values;  and  (c)  schematics  qualitatively  showing  the  evolution  of  the  film  thickness  as  well  as  the 
surface  roughness  of  the  CNT-AI2O3  tubes. 
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Young's  Modulus  (GPa) 


Figure  S2.  Measured  Young’s  modulus  (in  compression)  as  versus  micropillar  diameter  for  bare 
CNTs,  and  CNTs  coated  with  13,  30,  and  51  nm  AI2O3  by  ALD. 


24 


Compressive  modulus  [GPa] 


Figure  S3.  Plot  showing  the  sensitivity  of  the  model,  used  to  predict  Young’s  modulus,  to  the 
penetration  depth  of  the  ALD  coating  (hence  changing  the  coated  volume  of  the  micropillar). 
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Figure  S4.  Load-displacement  curves  for  AI2O3/CNT  micropillars,  with  number  of  ALD  cycles 
as  noted.  The  circle  superimposed  on  each  curve  indicates  the  strength  value  used  to  construct 
the  strength-density  relationship  in  Figure  6. 
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